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1  Statement  of  the  Problem  Studied 


The  objective  of  the  research  was  to  obtain  a  fundamental  understanding  of  the  physical 
and  chemical  mechanisms  of  autoignition  and  combustion  of  diesel  and  JP-8  in  non  premixed 
systems.  Experimental,  numerical  and  analytical  studies  were  carried  out.  Diesel  and  JP-8 
comprise  hundreds  of  aliphatic  and  aromatic  hydrocarbon  compounds.  The  major  components 
of  these  fuels  are  straight  chain  paraffins,  branched  chain  paraffins,  cycloparaffins,  aromatics, 
and  alkenes  [1-9].  A  composition  of  diesel  is  given  in  Refs.  [4]  and  [6]  showing  by  volume 
roughly  41  %  paraffins,  30  %  cycloparaffins,  and  29  %  aromatics.  The  composition  of  diesel 
given  in  weight  percent  in  Ref.  [3]  is  39%  paraffins,  44%  cycloparaffins,  and  17%  aromatics. 
In  JP-8  the  concentration  of  paraffins  is  on  the  average  60  %  by  volume,  that  of  cycloparaffins 
20%,  that  of  aromatics  18%,  and  that  alkenes  2%  [1].  Detailed  chemical  kinetic  mechanisms 
describing  combustion  for  many  of  the  components  in  diesel  as  well  as  those  in  JP-8  are  not 
available.  Therefore,  an  objective  of  this  research  was  to  develop  surrogates  for  JP-8  and 
diesel.  These  surrogates  are  mixtures  of  several  hydrocarbon  compounds. 

Edwards  and  Maurice  [1]  have  defined  two  types  of  surrogates:  physical  surrogate  and 
chemical  surrogate.  Physical  surrogates  are  mixtures  that  have  generally  the  same  physical 
properties  as  diesel  or  JP-8.  Physical  surrogates  are  unlikely  to  reproduce  aspects  of  com¬ 
bustion  of  diesel  or  JP-8.  Chemical  surrogates  are  mixtures  that  have  generally  the  same 
chemical-class  composition  and  average  molecular  weight  as  diesel  or  JP-8  [1].  In  general 
chemical  surrogates  are  designed  to  reproduce  aspects  of  combustion  of  diesel  or  JP-8.  The 
type  of  surrogate  would  depend  on  the  characteristics  being  simulated.  Edwards  and  Mau¬ 
rice  [1]  have  suggested  guidelines  for  developing  surrogates.  The  objective  of  this  study  was 
to  develop  chemical  surrogates  that  accurately  reproduce  autoignition  and  extinction  charac¬ 
teristics  of  diesel  and  JP-8  in  nonpremixed  systems. 

Combustion  processes  in  diesel  engines  closely  resemble  non  premixed  systems.  Therefore, 
the  research  was  focused  on  non  premixed  flames.  In  non  premixed  systems  the  flow  field 
and  combustion  chemistry  of  the  fuels  have  a  significant  inffuence  on  the  critical  conditions  of 
autoignition  and  flame  extinction.  Their  influence  can  be  discussed  in  terms  of  the  Damkohler 
number,  6,  defined  as  the  ratio  of  a  characteristic  ffow  time,  t/,  to  a  characteristic  chemical 
reaction  time,  rc.  The  reciprocal  of  the  flow  time  can  be  described  either  by  the  strain  rate, 
a,  or  the  scalar  dissipation  rate,  x  [10-12],  The  characteristic  chemical  time  depends  on  the 
type  of  fuel  and  the  initial  temperature  of  the  reactants.  Previous  analytical  studies  on  non 
premixed  systems  have  shown  that  autoignition  will  take  place  if  the  value  of  6  is  greater  than 
some  critical  value  di  and  flame  extinction  will  take  place  if  the  value  of  6  is  less  than  some 
critical  value  5e  [10-14],  Therefore  for  fixed  rc,  autoignition  will  take  place  for  ry  greater  than 
some  critical  value  or  for  the  strain  rate  and  the  scalar  dissipation  rate  less  than  some  critical 
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value,  represented  by  ai  and  Xh  respectively.  For  fixed  rc,  the  flame  will  extinguish  for  Tf 
less  than  some  critical  value  or  for  the  strain  rate  and  the  scalar  dissipation  rate  greater  than 
some  critical  value,  represented  by  ae  and  Xe>  respectively. 

The  research  was  performed  using  paraffins  (r?.-heptane,  n-decane,  and  n-dodecane)  cy¬ 
cloparaffins  (methylcyclohexane,  cyclohexane),  and  aromatics  (toluene,  o-xylene)  because  they 
represent  the  types  of  fuels  in  diesel  and  JP-8.  Detailed  chemical  kinetic  mechanisms  describ¬ 
ing  oxidation  of  a  number  of  these  compounds  are  either  already  available  or  they  are  likely 
to  be  developed.  Experimental  data  on  autoignition  and  flame  extinction  were  obtained.  The 
experimental  data  was  used  to  validate  detailed  chemical  kinetic  mechanisms. 

The  research  results  are  described  in  the  journal  articles  1-4,  listed  in  section  3.1,  and  in 
the  technical  reports  1 — 10  listed  in  section  3.2.  A  number  of  these  technical  reports  will  be 
revised  and  submitted  for  publications  in  journals.  Summaries  of  the  results  of  the  research 
are  described  in  section  2. 

2  Summary  of  the  Most  Important  Results 

The  summaries  of  the  important  results  are  given  in  this  section.  The  research  described  in 
sections  2.1,  2.2,  2.3,  2.5,  and  2.6  are  conducted  on  pure  fuels.  The  results  are  useful  for  de¬ 
veloping  surrogate  fuels.  It  has  been  established  that  that  the  chemical-kinetic  mechanisms  of 
oxidation  of  low  molecular  weight  fuels  such  as  hydrogen,  carbon  monoxide,  methane,  ethene 
(C2H4),  ethane  (C2Hq),  propene  (C3H6),  and  propane  (CsHg)  are  subsets  of  chemical-kinetic 
mechanisms  of  oxidation  of  n-heptane,  n-decane,  and  n-hexadecane.  Therefore  accurate  de¬ 
scriptions  of  mechanisms  of  oxidation  of  these  low  molecular  weight  fuels  are  necessary  first 
steps  in  the  development  of  chemical-kinetic  mechanisms  of  oxidation  of  practical  fuels.  Accu¬ 
rate  experimental  data  obtained  over  a  wide  range  of  conditions  and  different  configurations 
are  required  to  test  the  accuracies  of  chemical-kinetic  mechanisms  of  combustion  of  various 
fuels.  Therefore  section  2.4  gives  experimental  data  that  can  be  used  to  test  predictions  of  ex¬ 
tinction  and  autoignition  of  chemical-kinetic  mechanisms  describing  oxidation  of  C2H4,  C2H6, 
C3H6,  and  C3H8.  Section  2.7  describes  studies  leading  to  the  development  of  surrogate  fuels. 
Three  body  reactions,  for  example  H  +  O2  +  M  — >  HO2  +  M,  where  M  represents  any  third 
body  are  key  steps  in  chemical-kinetic  mechanisms  of  oxidation  of  hydrocarbon  fuels.  The 
chaperon  efficiency  of  the  third  body  is  not  well  known.  An  experimental  and  numerical  study 
focused  on  the  influence  of  this  reaction  is  described  in  section  2.8. 
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2.1  Ignition  in  the  Viscous  Layer  Between  Counterflowing  Streams:  Asymp¬ 
totic  Theory  with  Comparison  to  Experiments 

This  research  is  described  in  journal  article  1  listed  in  section  3.1  and  the  report  1  listed  in 
section  3.2.  The  research  was  performed  in  collaboration  with  Professor  A.  Linan  at  Departa- 
rnento  de  Motopropulsion  y  Termofluidodinamica,  Universidad  Politecnica  de  Madrid,  Spain. 

A  formulation  is  given  for  describing  autoignition  in  non  premixed  systems.  Steady  lami¬ 
nar  flow  of  two  counterflowing  streams  toward  a  stagnation  plane  is  considered.  One  stream 
comprises  fuel  and  the  other  oxygen.  The  characteristic  Reynolds  numbers  of  the  counter¬ 
flowing  streams  are  presumed  to  be  large  so  that  the  thickness  of  the  viscous  layer  formed 
in  the  vicinity  of  the  stagnation  plane  is  small.  The  chemical  reaction  that  takes  place  in 
the  viscous  layer  between  fuel  and  oxygen  is  described  by  a  one-step  overall  process.  The 
activation  energy  of  the  reaction  is  presumed  to  be  large  in  comparison  to  the  thermal  en¬ 
ergy.  The  asymptotic  theory  developed  here  makes  available  explicit  formulas  for  predicting 
autoignition  in  the  viscous  layer.  From  these  results  a  simple  but  reasonably  accurate  method 
is  developed  for  deducing  the  activation  energy,  E,  and  frequency  factor,  B,  of  the  rate  of  the 
one-step  reaction  between  the  fuel  and  oxygen.  To  illustrate  the  application  of  this  method, 
experiments  are  carried  out  in  the  counterflow  configuration.  The  fuels  tested  are  n-heptane, 
n-decane,  JP-10,  and  toluene.  Experimental  data  obtained  are  the  velocities  and  tempera¬ 
tures  of  counterflowing  streams  at  autoignition.  Values  of  E  and  B  are  obtained  by  using  the 
experimental  data  in  the  formulas  given  by  the  asymptotic  theory.  These  values  of  E  and  B 
are  found  to  agree  well  with  those  obtained  from  numerical  calculations. 

2.2  Extinction  and  Autoignition  of  n-Heptane  in  Counterflow  Configuration 

This  research  is  described  in  journal  article  2  listed  in  section  3.1  and  the  report  2  listed  in 
section  3.2.  The  research  was  performed  in  collaboration  Dr.  W.  J.  Pitz  and  Dr.  H.  J.  Curran 
at  Lawrence  Livermore  National  Laboratory,  Livermore,  California. 

A  study  is  performed  to  elucidate  the  mechanisms  of  extinction  and  autoignition  of  n-heptane 
in  strained  laminar  flows  under  non  premixed  conditions.  A  previously  developed  detailed 
mechanism  made  up  of  2540  reversible  elementary  reactions  among  557  species  is  the  start¬ 
ing  point  for  the  study.  The  detailed  mechanism  was  previously  used  to  calculate  ignition 
delay  times  in  homogeneous  reactors,  and  concentration  histories  of  a  number  of  species  in 
plug-flow  and  jet-stirred  reactors.  An  intermediate  mechanism  made  up  of  1282  reversible 
elementary  reactions  among  282  species  and  a  short  mechanism  made  up  of  770  reversible 
elementary  reactions  among  160  species  are  assembled  from  this  detailed  mechanism.  Ignition 
delay  times  in  an  isochoric  homogeneous  reactor  calculated  using  the  intermediate  and  the 
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short  mechanism  are  found  to  agree  well  with  those  calculated  using  the  detailed  mechanism. 
The  intermediate  and  the  short  mechanism  are  used  to  calculate  extinction  and  autoignition 
of  n-heptane  in  strained  laminar  flows.  Steady  laminar  flow  of  two  counterflowing  streams 
toward  a  stagnation  plane  is  considered.  One  stream  made  up  of  prevaporized  n-heptane  and 
nitrogen  is  injected  from  the  fuel  boundary  and  the  other  stream  made  up  of  air  and  nitrogen 
is  injected  from  the  oxidizer  boundary.  Critical  conditions  of  extinction  and  autoignition  given 
by  the  strain  rate,  temperature  and  concentrations  of  the  reactants  at  the  boundaries,  are  cal¬ 
culated.  The  results  are  found  to  agree  well  with  experiments.  Sensitivity  analysis  is  carried 
out  to  evaluate  the  influence  of  various  elementary  reactions  on  autoignition.  At  all  values 
of  the  strain  rate  investigated  here,  high  temperature  chemical  processes  are  found  to  control 
autoignition.  In  general,  the  influence  of  low  temperature  chemistry  is  found  to  increase  with 
decreasing  strain.  A  key  finding  of  the  present  study  is  that  strain  has  more  influence  on  low 
temperature  chemistry  than  the  temperature  of  the  reactants. 

2.3  Temperature  Cross-Over  and  Non-Thermal  Runaway  at  Two-Stage  Ig¬ 
nition  of  N-Heptane 

This  research  is  described  in  journal  article  3  listed  in  section  3.1  and  the  report  4  listed  in 
section  3.2.  The  research  was  performed  in  collaboration  with  Professor  N.  Peters  at  Institut 
fur  Technische  Mechanik,  RWTH  Aachen,  52056  Aachen,  Germany. 

To  calculate  ignition  delay  times  a  skeletal  56-step  mechanism  for  n-heptane  is  further 
reduced  to  a  short  30-step  mechanism  containing  two  isomers  of  the  n-heptyl-radical  and 
reactions  describing  both  the  high  temperature  and  the  low  temperature  chemistry.  This 
mechanism  reproduces  ignition  delay  times  at  various  pressures  and  temperatures  reasonably 
well.  Steady  state  assumptions  for  many  of  the  intermediate  species  are  introduced  to  derive 
separately  two  global  mechanisms  for  the  low  temperature  regime  as  well  as  for  the  interme¬ 
diate  and  high  temperature  regime.  In  those  formulations  the  OH  radical  is  depleted  by  fast 
reactions  with  the  fuel,  as  long  as  fuel  is  present.  Its  steady  state  relation  shows  that  the  OH 
concentration  would  blow  up  as  soon  as  the  fuel  is  depleted.  Therefore  the  depletion  of  the 
fuel  is  used  as  a  suitable  criterion  for  ignition. 

In  the  intermediate  temperature  regime  the  first  stage  ignition  is  related  to  a  change  from 
chain-branching  to  chain-breaking  as  the  temperature  crosses  a  certain  threshold.  The  chain 
branching  reactions  result  in  a  build-up  of  ketohydroperoxides  which  dissociate  to  produce 
OH  radicals.  This  is  associated  with  a  slight  temperature  rise  which  leads  to  a  crossing  of 
the  threshold  temperature  with  the  consequence  that  the  production  of  OH  radicals  by  keto¬ 
hydroperoxides  suddenly  ceases.  The  subsequent  second  stage  is  driven  by  the  much  slower 
production  of  OH  radicals  owing  to  the  dissociation  of  hydrogen  peroxide.  The  OH  radicals 
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react  with  the  fuel  at  nearly  constant  temperature  until  the  latter  is  fully  depleted. 

In  all  three  regimes  analytical  solutions  for  the  ignition  delay  time  are  presented.  The  re¬ 
duced  4-step  mechanism  of  the  low  temperature  regime  leads  with  the  assumption  of  constant 
temperature  to  linear  differential  equations,  which  are  solved.  The  calculated  ignition  delay 
times  at  fuel  depletion  compares  well  with  those  of  the  30-step  mechanism.  The  analysis  for 
the  intermediate  temperature  regime  starts  from  a  4-step  subset  of  a  9-step  reduced  mecha¬ 
nism.  It  contains  the  cross-over  dynamics  in  form  of  a  temperature  dependent  stoichiometric 
coefficient  which  is  analyzed  mathematically.  The  resulting  closed  form  solutions  describe  the 
first  stage  ignition,  the  temperature  cross-over  and  the  second  stage  ignition.  They  also  iden¬ 
tify  the  rate  determining  reactions  and  quantify  the  influence  of  their  rates  on  the  first  and  the 
second  ignition  times.  The  high  temperature  regime  is  governed  by  a  three-step  mechanism 
leading  to  a  nonlinear  problem  which  is  solved  by  asymptotic  analysis. 

While  the  dissociation  reaction  of  the  ketohydroperoxide  dominates  the  low  temperature 
regime  and  the  first  stage  ignition  of  the  intermediate  temperature  regime,  the  hydrogen 
peroxide  dissociation  takes  this  role  for  the  second  stage  of  the  intermediate  and  in  the  high 
temperature  regime.  The  overall  activation  energy  of  the  ignition  delay  time  in  the  low 
temperature  regime  is  the  mean  of  the  activation  energies  of  two  reactions  only.  The  overall 
activation  energy  of  the  ignition  delay  time  in  the  high  temperature  regime  is  shown  to  be 
related  to  the  activation  energies  of  only  three  but  different  rate  determining  reactions. 

2.4  Nonpremixed  and  Premixed  Extinction  and  Autoignition  of  C2H1,  C2H6, 
C3H6,  and  C3H8 

This  research  is  described  in  journal  article  4  listed  in  section  3.1  and  report  3  and  report  6 
listed  in  section  3.2. 

Experimental  studies  are  conducted  on  laminar  nonpremixed  and  premixed  flames  stabi¬ 
lized  in  the  counterflow  configuration.  The  fuels  tested  are  ethene  (C2H4),  ethane  (C2H6), 
propene  (C3Hq),  and  propane  (C3Hg).  Studies  on  nonpremixed  systems  are  carried  out  by 
injecting  a  fuel  stream  made  up  of  fuel  and  nitrogen  (N2)  from  one  duct  and  an  oxidizer 
stream  made  up  of  air  and  N2  from  the  other  duct.  Studies  on  premixed  systems  are  carried 
out  by  injecting  a  premixed  reactant  stream  made  up  of  fuel,  oxygen,  and  nitrogen  from  one 
duct,  and  an  inert-gas  stream  of  N2  from  the  other  duct.  Critical  conditions  of  extinction 
are  measured  by  increasing  the  flow  rates  of  the  counterflowing  streams  until  the  flame  extin¬ 
guishes.  Critical  conditions  of  autoignition  are  measured  by  preheating  the  oxidizer  stream 
of  the  non  premixed  system  and  the  inert-gas  stream  of  the  premixed  system.  Experimental 
data  for  autoignition  are  obtained  over  a  wide  range  of  temperatures  of  the  heated  streams.  In 
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addition  for  premixed  systems  experimental  data  are  obtained  for  a  wide  range  of  values  of  the 
equivalence  ratio  including  fuel-lean  and  fuel-rich  conditions.  Numerical  calculations  are  per¬ 
formed  using  a  detailed  chemical-kinetic  mechanism  and  compared  with  measurements.  The 
present  study  highlights  the  influences  of  nonuniform  flow  field  on  autoignition  in  premixed 
systems  that  were  not  available  from  previous  studies  using  shock  tubes.  For  the  premixed 
system  considered  here  the  changes  in  the  strain  rates  at  extinction  with  equivalence  ratio 
are  found  to  be  similar  to  previous  observations  of  changes  in  laminar  burning  velocities  with 
equivalence  ratio.  The  studies  on  autoignition  in  the  premixed  system  show  that  the  tem¬ 
perature  of  the  inert-gas  stream  at  autoignition  reaches  a  minimum  for  a  certain  equivalence 
ratio.  For  premixed  systems  abrupt  extinction  and  autoignition  are  not  observed  if  the  value 
of  the  equivalence  ratio  is  less  than  some  critical  value. 

2.5  Chemical  Kinetic  Study  of  Toluene  Oxidation 

This  research  is  described  in  report  5  and  report  7  listed  in  section  3.2.  The  research  was 
performed  in  collaboration  with  W.  J.  Pitz  and  C.  K.  Westbrook  at  Lawrence  Livermore  Na¬ 
tional  Laboratory,  Livermore,  California,  J.  W.  Bozzelli,  and  C.-J.  Chen,  at  Chemistry  and 
Chemical  Engineering  Department,  New  Jersey  Institute  of  Technology,  Newark,  NJ,  and  I. 
Da  Costa,  R.  Fournet,  F.  Billaud,  F.  Battin-Leclerc,  at  Departement  de  Chime  Physique  des 
Reactions,  CNRS-ENSIC,  BP.  451,  1,  rue  Grandville,  51001  Nancy,  France. 

A  study  was  performed  to  elucidate  the  chemical-kinetic  mechanism  of  combustion  of 
toluene.  A  detailed  chemical-kinetic  mechanism  for  toluene  was  improved  by  adding  a  more 
accurate  description  of  the  phenyl  +  O2  reaction  channels,  toluene  decomposition  reactions 
and  the  benzyl  +  O  reaction.  Results  of  the  chemical  kinetic  mechanism  are  compared  with 
experimental  data  obtained  from  premixed  and  non  premixed  systems.  Under  premixed  con¬ 
ditions,  predicted  ignition  delay  times  are  compared  with  new  experimental  data  obtained 
in  shock  tube.  Also,  calculated  species  concentration  histories  are  compared  to  experimental 
flow  reactor  data  from  the  literature.  Under  non  premixed  conditions,  critical  conditions  of 
extinction  and  autoignition  were  measured  in  strained  laminar  flows  in  the  counterflow  con¬ 
figuration.  Numerical  calculations  are  performed  using  the  chemical-kinetic  mechanism  at 
conditions  corresponding  to  those  in  the  experiments.  Critical  conditions  of  extinction  and 
autoignition  are  predicted  and  compared  with  the  experimental  data.  Comparisons  between 
the  model  predictions  and  experimental  results  of  ignition  delay  times  in  shock  tube,  and  ex¬ 
tinction  and  autoignition  in  non  premixed  systems  show  that  the  chemical-kinetic  mechanism 
predicts  that  toluene/air  is  overall  less  reactive  than  observed  in  the  experiments.  For  both 
premixed  and  non  premixed  systems,  sensitivity  analysis  was  used  to  identify  the  reaction 
rate  constants  that  control  the  overall  rate  of  oxidation  in  each  of  the  systems  considered. 
Under  shock  tube  conditions,  the  reactions  that  influence  ignition  delay  time  are  H  +  O2 


chain  branching,  the  toluene  decomposition  reaction  to  give  an  H  atom,  and  the  toluene  + 
H  abstraction  reaction.  The  reactions  that  influence  autoignition  in  non  premixed  systems 
involve  the  benzyl  +  HO2  reaction  and  the  phenyl  +  O2  reaction. 

Further  details  of  this  work  is  given  in  Appendix  A.  A  journal  publication  describing 
results  of  this  research  is  under  preparation. 

2.6  The  Structure  of  Nonpremixed  n-Decane  Flames 

This  research  is  described  in  report  8  listed  in  section  3.2. 

An  experimental  and  numerical  study  is  performed  to  elucidate  the  structure  and  mech¬ 
anisms  of  extinction  and  autoignition  n-decane  flames  in  strained  laminar  flows  under  non¬ 
premixed  conditions.  Experiments  are  conducted  on  flames  stabilized  between  two  counter¬ 
flowing  streams.  The  fuel-stream  is  a  mixture  of  prevaporized  n-decane  and  nitrogen,  and 
the  oxidizer  stream  is  a  mixture  of  air  and  nitrogen.  Concentration  profiles  of  C10H22,  O2, 
N2,  CO2,  CO,  H2O,  CO,  CH4,  C2-hydrocarbons,  C3-hydrocarbons,  and  C4-hydrocarbons,  are 
measured.  The  measurements  are  made  by  removing  gas  samples  from  the  flame  using  a  quartz 
microprobe  and  analyzing  the  samples  using  a  gas  chromatograph.  Temperature  profiles  are 
measured  using  a  thermocouple.  Numerical  calculations  are  performed  using  detailed  chem¬ 
istry  to  determine  the  flame  structure  and  the  results  are  compared  with  the  measurements. 
Critical  conditions  of  extinction  are  measured.  Data  giving  the  mass  fraction  of  oxygen  in 
the  oxidizer  stream  as  a  function  of  the  strain  rate  at  extinction  are  obtained.  Critical  condi¬ 
tions  of  autoignition  are  measured.  Data  giving  the  temperature  of  the  oxidizer  stream  as  a 
function  of  the  strain  rate  are  obtained.  Numerical  calculations  are  performed  using  detailed 
chemistry.  Critical  conditions  of  extinction  and  autoignition  are  obtained  and  compared  with 
measurements. 

Further  details  of  this  work  is  given  in  Appendix  B.  A  journal  publication  describing 
results  of  this  research  is  under  preparation. 

2.7  Autoignition  and  Extinction  of  Hydrocarbon  Fuels  in  Nonpremixed 
Systems 

This  research  is  described  in  report  9  listed  in  section  3.2. 

Experimental  studies  are  carried  out  on  extinction  and  autoignition  of  liquid  hydrocarbon 
fuels  in  nonpremixed  systems.  The  counterflow  configuration  is  employed.  In  this  configura¬ 
tion  an  oxidizer  made  up  of  air  and  nitrogen  flows  over  the  surface  of  a  pool  of  liquid  fuel. 
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Critical  conditions  of  extinction  are  measured  for  n-heptane,  iso-octane,  n-octane,  n-decane, 
n-dodecane,  n-hexadecane,  methylcyclohexane,  cyclohexane,  diesel  and  JP-8.  Data  giving  the 
mass  fraction  of  oxygen  in  the  oxidizer  stream  as  a  function  of  the  strain  rate  at  extinction 
are  obtained.  Critical  conditions  of  autoignition  are  measured  for  n-heptane,  iso-octane,  n- 
octane,  n-decane,  n-dodecane,  n-hexadecane,  methylcyclohexane,  o-xylene,  JP-10,  diesel  and 
JP-8.  The  experimental  results  show  that  the  critical  conditions  of  extinction  of  iso-octane, 
hexadecane,  and  a  mixture  of  70%  decane  and  30%  xylene  are  similar  to  that  of  diesel/  The 
critical  conditions  of  autoignition,  however,  of  these  fuels  and  fuel  mixture  are  not  similar  to 
that  of  diesel.  This  illustrates  the  challenges  in  developing  surrogate  fuels. 

Further  details  of  this  work  is  given  in  Appendix  C.  A  journal  publication  describing 
results  of  this  research  is  under  preparation. 

2.8  The  Influence  of  Water  Vapor  of  the  Critical  Conditions  of  Extinction 
and  Autoignition  of  Premixed  Ethene  Flames 

This  research  is  described  in  report  10  listed  in  section  3.2. 

An  experimental  and  numerical  study  is  performed  to  elucidate  the  influence  of  addition  of 
water  vapor  on  the  structure  and  mechanisms  of  extinction  and  autoignition  premixed  ethene 
(C2H4)  flames  in  strained  laminar  flows.  Experiments  are  conducted  on  flames  stabilized 
between  two  counterflowing  streams.  The  counterflow  burner  is  made  up  of  two  ducts.  A 
premixed  reactant  stream  made  up  of  C2H4,  air,  and  nitrogen  is  injected  from  one  duct,  and 
an  inert  gas  stream  of  N2  is  injected  from  the  other  duct.  Water  vapor  is  added  to  the  pre¬ 
mixed  reactant  stream  using  a  vaporizer.  The  vaporizer  contains  liquid  water  that  is  heated 
using  an  electrical  heater.  The  temperature  of  water  is  measured  using  a  thermocouple.  The 
water- vapor  in  the  vaporizer  is  presumed  to  be  in  equilibrium  with  the  liquid.  A  part  of  the 
air  in  premixed  reactant  stream  percolates  through  the  liquid  water  in  the  vaporizer.  The 
mass  fraction  of  water  vapor  in  the  premixed  reactant  stream,  1h2Oj  can  be  changed  either 
by  changing  the  temperature  of  the  water  in  the  vaporizer  or  the  flow  rate  of  air  through 
the  vaporizer.  The  premixed  reactant  stream  is  characterized  by  the  equivalence  ratio  cj)  and 
the  adiabatic  temperature  Ta(\.  To  clarify  the  influence  of  water  vapor  on  extinction  and 
autoignition  studies  are  carried  out  at  fixed  values  of  (j)  and  Tad.  Critical  conditions  of  extinc¬ 
tion  are  measured  by  increasing  the  flow  rates  of  the  counterflowing  streams  until  the  flame 
extinguishes.  Critical  conditions  of  autoignition  are  measured  by  preheating  the  inert  gas 
stream  of  the  premixed  system.  Experimental  data  are  obtained  for  a  wide  range  of  values  of 
Yh20-  Numerical  calculations  are  performed  using  a  detailed  chemical  kinetic  mechanism  and 
compared  with  measurements. 
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Further  details  of  this  work  is  given  in  Appendix  D.  A  journal  publication  describing 

results  of  this  research  is  under  preparation. 

3  List  of  all  Publications  and  Technical  Reports 

3.1  Journal  Articles 

1.  Seiser,  R.,  Seshadri,  K.,  Piskernik,  E.,  and  Lihan,  A.,  “Ignition  in  the  Viscous  Layer  Be¬ 
tween  Counterflowing  Streams:  Asymptotic  Theory  with  Comparisons  to  Experiments” 
Combustion  and  Flame,  122,  2000,  pp  339-349. 

2.  Seiser,  R.,  Pitsch,  H.,  Seshadri,  K.,  Pitz,  W.  J.,  Curran,  H.  J.,  “Extinction  and  Au¬ 
toignition  of  n-Heptane  in  Counterflow  Configuration,”  Proceedings  of  the  Combustion 
Institute,  28,  2000,  pp  2029-2037. 

3.  Peters,  N.,  Paczko,  G.,  Seiser,  R.,  and  Seshadri,  K.,  “Temperature  Cross-Over  and  Non- 
Thermal  Runway  at  the  Two-Stage  Ignition  of  N-Heptane,”  Combustion  and  Flame, 
128,  2002,  pp  38-59. 

4.  Hurner,  S.,  Seiser,  R.,  and  Seshadri,  K.,  “Nonpremixed  and  Premixed  Extinction  and 
Autoignition  of  C2H4,  C2H6,  C3H6,  and  C3H8,”  Proceedings  of  the  Combustion  Institute, 
29,  2002,  accepted  for  publication. 

3.2  Papers  Published  in  Proceedings  of  Technical  Meetings 

1.  R.  Seiser,  K.  Seshadri,  E.  Piskernik,  and  A.  Linan,  “Ignition  in  the  Viscous  layer  Between 
Counterflowing  Streams:  Asymptotic  Theory  with  Comparison  to  Experiments,”  Paper 
ff  00S-27,  Western  States  Section  of  the  Combustion  Institute,  Spring  2000  Meeting, 
Colorado  School  of  Mines,  Golden,  CO  80401,  March  13 — 14,  2000. 

2.  R.  Seiser,  H.  Pitsch,  K.  Seshadri,  W.  J.  Pitz,  and  H.  J.  Curran,  “Extinction  and  Au¬ 
toignition  of  n-Heptane  in  Counterflow  Configuration,”  Paper  ff  00S-28,  Western  States 
Section  of  the  Combustion  Institute,  Spring  2000  Meeting,  Colorado  School  of  Mines, 
Golden,  CO  80401,  March  13—14,  2000. 

3.  R.  Seiser,  K.  Seshadri,  and  W.  J.  Pitz,  “Experimental  and  Numerical  Studies  of  Ex¬ 
tinction  and  Autoignition  of  C3H8,  C3H6,  C2H6,  and  C2H4,”  paper  ff  149,  2nd  Joint 
Meeting  of  the  US  Sections  of  the  Combustion  Institute,  Oakland  California,  March 
25—28,  2001. 

4.  N.  Peters,  G.  Paczko,  R.  Seiser,  and  K.  Seshadri,  “Temperature  Cross-Over  and  Non- 
Thermal  Runaway  at  Two-Stage  Ignition  of  n-Heptane,”  paper  ff  146,  2nd  Joint  Meeting 
of  the  US  Sections  of  the  Combustion  Institute,  Oakland  California,  March  25 — 28,  2001. 
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5.  W.  J.  Pitz  ,  C.  K.  Westbrook,  R.  Seiser,  K.  Seshadri,  J.  W.  Bozzelli,  I.  Da  Costa, 
R.  Fournet,  F.  Billaud,  and  F.  Battin-Leclerc,  “Chemical-Kinetic  Characterization  of 
Combustion  of  Toluene,”  “Chemical-Kinetic  Study  of  Toluene  Oxidation,”  paper  #  252, 
2nd  Joint  Meeting  of  the  US  Sections  of  the  Combustion  Institute,  Oakland  California, 
March  25—28,  2001. 

6.  S.  Hurner  ,  R.  Seiser,  and  K.  Seshadri,  “Counterflow  Ignition  of  Premixed  Ethene/Air 
Mixtures  with  Heated  Nitrogen,”  paper  #  01F-6,  Western  States  Section  of  the  Com¬ 
bustion  Institute,  Fall  2001  Meeting,  University  of  Utah,  Salt  Lake  City,  Utah,  October 
15—16,  2001. 

7.  W.  J.  Pitz,  R.  Seiser,  J.  W.  Bozzelli,  I.  Da  Costa,  R.  Fournet,  F.  Billaud,  F.  Battin- 
Leclerc,  K.  Seshadri,  and  C.  K.  Westbrook,  “Chemical-Kinetic  Study  of  Toluene  Oxida¬ 
tion,”  paper  #  OIF-28,  Western  States  Section  of  the  Combustion  Institute,  Fall  2001 
Meeting,  University  of  Utah,  Salt  Lake  City,  Utah,  October  15 — 16,  2001. 

8.  Tanoue,  K.,  Seiser,  R.,  and  Seshadri,  K.,  “The  Structure  of  Nonpremixed  n- Decane 
Flames,”  Paper  A31,  3rd  Joint  Meeting  of  the  United  States  Sections  of  the  Combustion 
Institute,  The  University  of  Illinois  at  Chicago,  Chicago,  March  16-19,  2003. 

9.  Hurner,  S.,  Seiser,  R.  ,  and  Seshadri,  K.,  “Autoignition  and  Extinction  of  Hydrocarbon 
Fuels  in  Nonpremixed  Systems,”  Paper  B26,  3rd  Joint  Meeting  of  the  United  States 
Sections  of  the  Combustion  Institute,  The  University  of  Illinois  at  Chicago,  Chicago, 
March  16-19,  2003. 

10.  Geieregger,  M.,  Seiser,  R.,  and  Seshadri,  K.,  “The  Influence  of  Water  Vapor  on  the 
Critical  Conditions  of  Extinction  and  Autoignition  pf  Premixed  Ethene  Flames,”  Paper 
B27,  3rd  Joint  Meeting  of  the  United  States  Sections  of  the  Combustion  Institute,  The 
University  of  Illinois  at  Chicago,  Chicago,  March  16-19,  2003. 


4  Participating  Scientific  Personnel 

The  participating  scientific  personnel  were  1)  Professor  K.  Seshadri,  2)  Dr.  R.  Seiser  (research 
staff),  3)  Dr.  K.  Tanoue  (visiting  scholar),  4)  Mr.  S.  Hurner  (graduate  student),  5)  Mr.  E. 
Piskernick  (graduate  student),  and  6)  M.  Geieregger  (graduate  student). 

5  Report  of  Inventions 

None. 
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Abstract 

A  study  was  performed  to  elucidate  the  chemical-kinetic  mechanism  of  combustion  of 
toluene.  A  detailed  chemical-kinetic  mechanism  for  toluene  was  improved  by  adding  a 
more  accurate  description  of  the  phenyl  +  O2  reaction  channels,  toluene  decomposition 
reactions  and  the  benzyl  +  O  reaction.  Results  of  the  chemical  kinetic  mechanism  are  com¬ 
pared  with  experimental  data  obtained  from  premixed  and  nonpremixed  systems.  Under 
premixed  conditions,  predicted  ignition  delay  times  are  compared  with  new  experimental 
data  obtained  in  shock  tube.  Also,  calculated  species  concentration  histories  are  compared 
to  experimental  flow  reactor  data  from  the  literature.  Under  nonpremixed  conditions,  crit¬ 
ical  conditions  of  extinction  and  autoignition  were  measured  in  strained  laminar  flows  in 
the  counterflow  configuration.  Numerical  calculations  are  performed  using  the  chemical- 
kinetic  mechanism  at  conditions  corresponding  to  those  in  the  experiments.  Critical  con¬ 
ditions  of  extinction  and  autoignition  are  predicted  and  compared  with  the  experimental 
data.  Comparisons  between  the  model  predictions  and  experimental  results  of  ignition 
delay  times  in  shock  tube,  and  extinction  and  autoignition  in  nonpremixed  systems  show 
that  the  chemical-kinetic  mechanism  predicts  that  toluene/air  is  overall  less  reactive  than 
observed  in  the  experiments.  For  both  premixed  and  nonpremixed  systems,  sensitivity 
analysis  was  used  to  identify  the  reaction  rate  constants  that  control  the  overall  rate  of 
oxidation  in  each  of  the  systems  considered.  Under  shock  tube  conditions,  the  reactions 
that  influence  ignition  delay  time  are  H  +  O2  chain  branching,  the  toluene  decomposition 
reaction  to  give  an  H  atom,  and  the  toluene  +  H  abstraction  reaction.  The  reactions  that 
influence  autoignition  in  nonpremixed  systems  involve  the  benzyl  +  HO2  reaction  and  the 
phenyl  +  O2  reaction. 


Introduction 

Alkylated  benzenes  are  an  important  class  of  hydrocarbons  because  they  comprise  a  signifi¬ 
cant  portion  of  gasoline  and  diesel  fuels.  Knowledge  of  the  oxidation  chemistry  of  alkylated 
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benzenes  is  needed  in  developing  predictive  models  that  can  treat  autoignition,  and  premixed 
and  nonpremixed  burning  of  transportation  fuels  in  internal  combustion  engines.  Toluene 
(C6H5CH3)  has  one  of  the  simplest  molecular  structures  of  the  alkylated  benzenes  and  is  a 
reasonable  starting  point  for  the  development  of  detailed  chemical-kinetic  reaction  mecha¬ 
nisms  for  alkylated  benzenes.  Much  previous  work  has  been  done  on  the  oxidation  of  toluene. 
Several  research  groups  have  developed  detailed  chemical-kinetic  reaction  mechanisms  for 
toluene.  Most  recently,  Klotz  et  al.  [1]  supplemented  the  toluene  mechanism  of  Erndee  et 
al.  [2]  to  improve  the  predictions  for  the  intermediates  1,3  butadiene,  acetylene  and  benzalde- 
hyde.  Zhong  and  Bozzelli  [3-5]  developed  a  more  accurate  description  of  radical  additions 
to  cyclopentadiene  and  associations  with  cyclopentadienyl  radical;  they  included  these  reac¬ 
tions  in  a  detailed  chemical-kinetic  mechanism  for  toluene  that  they  developed.  Lindstedt 
and  Maurice  [6]  developed  a  very  comprehensive  toluene  mechanism  whose  predictions  they 
compared  to  experimental  results  from  counterflow  diffusion  flames,  plug  flow  reactors,  shock 
tubes  and  premixed  flames.  Erndee  et  al.  [2]  developed  a  detailed  chemical- kinetic  mechanism 
for  toluene  that  was  benchmark  for  many  years. 

There  are  quite  a  few  experimental  studies  of  toluene  oxidation  whose  data  are  very  useful 
for  mechanism  validation.  Several  experimental  studies  of  toluene  oxidation  in  a  flow  reactor 
were  performed  at  Princeton  University  [1,  2,  7,  8].  Ignition  of  toluene  in  a  rapid  compression 
machine  was  performed  by  Griffiths  et  al.  [9]  and  by  Roubaud  et  al.  [10].  Their  rapid  com¬ 
pression  machine  results  show  that  toluene  oxidation  chemistry  lacks  the  two  stage  ignition 
observed  in  paraffinic  fuels.  Experimental  data  for  the  critical  conditions  of  autoignition  of 
toluene  in  the  counterflow  configuration  are  given  in  Ref.  [11].  Using  this  experimental  data, 
overall  chemical-kinetic  rate  parameters  that  characterize  the  rate  of  one-step  overall  reaction 
between  fuel  and  oxygen  were  obtained  [11].  The  present  work  offers  new  results  for  shock 
tube  ignition  of  toluene  and  new  results  for  extinction  under  nonpremixed  conditions.  These 
experiments  provide  additional  data  for  validating  chemical-kinetic  mechanisms  for  toluene. 
In  particular,  the  results  under  nonpremixed  conditions  highlight  the  influence  of  strain  on 
autoignition  and  extinction. 


Experimental 

Shock  tube  experiments 

The  experiments  were  performed  in  a  stainless  steel  78  mm  diameter  shock  tube  at  DCPR- 
CNRS-NANCY.  The  reaction  and  driver  parts  were  separated  by  two  terphane  diaphragms, 
which  were  ruptured  by  suddenly  decreasing  the  pressure  in  the  space  separating  them.  The 
driver  gas  was  helium  and  the  reacting  mixture  was  diluted  in  argon.  The  incident  and  re¬ 
flected  shock  velocities  were  measured  by  four  piezo-electric  pressure  transducers  located  along 
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the  reaction  section.  The  state  of  the  test  gas  behind  the  incident  and  the  reflected  shock  waves 
was  derived  from  the  value  of  the  incident  shock  velocity  by  using  classical,  one-dinrensional, 
shock  equations  of  mass,  momentum,  and  energy  conservation  applied  to  an  ideal  gas. 

The  pressure  profile  displayed  three  rises,  which  were  due  to  the  incident  shock  wave, 
the  reflected  shock  wave  and  the  ignition,  respectively.  The  onset  of  ignition  was,  however, 
most  accurately  detected  by  OH  radical  emission  at  306  nrn  through  a  quartz  window  with 
a  photomultiplier  fitted  with  a  monochromator  at  the  end  of  the  reaction  section.  The  last 
pressure  transducer  was  located  at  the  same  place  along  the  axis  of  the  tube  as  the  quartz 
window.  The  ignition  delay  time  was  defined  as  the  time  interval  between  the  pressure  rise 
measured  by  the  last  pressure  transducer  due  to  the  arrival  of  the  reflected  shock  wave  and  the 
rise  of  the  optical  signal  by  the  photomultiplier  up  to  10%  of  its  maximum  value.  The  igni¬ 
tion  delay  times  were  corrected  for  blast  wave  effects  by  adding  3  ps  to  the  measured  time  [12]. 

Oxygen  and  helium  were  purchased  from  Air  Liquide-Alphagaz  and  toluene  was  provided 
by  Aldrich  (purity:  99.8%).  The  toluene  concentration  in  the  reactants  was  kept  constant  at 
1.25  %.  The  values  of  equivalence  ratio,  (j),  examined  were  0.5,  1.0  and  1.5.  The  reflected  shock 
pressures  and  temperatures  ranged  from  8.0  to  9.4  atm  and  1300  to  1900  K,  respectively. 

Experiments  under  Nonpremixed  Conditions 

Experiments  under  nonpremixed  conditions  were  carried  out  in  the  counterflow  configuration. 
Figure  1  shows  a  schematic  illustration  of  the  counterflow  configuration.  Steady,  axisymmet- 
ric,  laminar  flow  of  two  counterflowing  streams  toward  a  stagnation  plane  is  considered.  In 
this  configuration  a  fuel  stream  made  up  of  prevaporized  toluene  and  nitrogen  is  injected  from 
the  fuel-duct,  and  an  oxidizer  stream  of  air  is  injected  from  the  oxidizer-duct.  These  jets 
flow  into  the  mixing  layer  between  the  two  ducts.  The  exit  of  the  fuel-duct  is  called  the  fuel 
boundary  and  the  exit  of  the  oxidizer-duct  the  oxidizer  boundary.  The  mass  fraction  of  fuel, 
the  temperature,  and  the  component  of  the  flow  velocity  normal  to  the  stagnation  plane  at 
the  fuel  boundary  are  represented  by  Yp,i,  Ti,  and  Vj ,  respectively.  The  mass  fraction  of  oxy¬ 
gen,  the  temperature,  and  the  component  of  the  flow  velocity  normal  to  the  stagnation  plane 
at  the  oxidizer  boundary  are  represented  by  Yo2p,  ^2,  and  V2,  respectively.  The  tangential 
components  of  the  flow  velocities  at  the  boundaries  are  presumed  to  be  equal  to  zero.  The 
distance  between  the  fuel  boundary  and  the  oxidizer  boundary  is  represented  by  L. 

In  the  experiments  the  momenta  of  the  counterflowing  reactant  streams  piVi2,i  =  1,2 
at  the  boundaries  are  kept  equal  to  each  other.  Here  p\  and  P2  represent  the  density  of 
the  mixture  at  the  fuel  boundary  and  at  the  oxidizer  boundary,  respectively.  This  condition 
ensures  that  the  stagnation  plane  formed  by  the  two  streams  is  approximately  in  the  middle 
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Figure  1:  Schematic  illustration  of  the  counterflow  configuration. 


of  the  region  between  the  two  boundaries.  The  value  of  the  strain  rate,  defined  as  the  normal 
gradient  of  the  normal  component  of  the  flow  velocity,  changes  from  the  fuel  boundary  to  the 
oxidizer  boundary  [13].  The  characteristic  strain  rate  on  the  oxidizer  side  of  the  stagnation 
plane  a  2  is  presumed  to  be  given  by  [13] 


a  2 


2\V2\(  \Vi\Vp~i\ 
L  \  |F2|  vW' 


(1) 


Equation  1  is  obtained  from  an  asymptotic  theory  where  the  Reynolds  numbers  of  the  lam¬ 
inar  flow  at  the  boundaries  are  presumed  to  be  large  [13].  Critical  conditions  of  extinction 
are  presumed  to  be  given  by  the  strain  rate,  <32, e)  and  the  mass  fraction  of  fuel  at  the  fuel 
boundary.  Critical  conditions  of  autoignition  are  presumed  to  be  given  by  the  strain  rate, 
02,7,  the  temperature  of  the  oxidizer  stream,  T2,/,  and  the  mass  fraction  of  fuel  and  at  the  fuel 
boundary.  The  experiments  were  conducted  at  a  pressure  of  1.013  bar. 


A  detailed  description  of  the  burner  is  given  elsewhere  [11, 14].  The  flow  rates  of  gases 
were  measured  by  computer-regulated  mass  flow  controllers.  The  velocities  of  the  reactants 
at  the  boundaries  were  presumed  to  be  equal  to  the  ratio  of  their  volumetric  flowrates  to  the 
cross-section  area  of  the  ducts.  The  temperature  of  the  fuel  stream  and  the  temperature  of 
the  oxidizer  stream  at  the  boundaries  were  measured  using  thermocouples.  The  measured 


18 


temperatures  were  corrected  for  radiative  heat  losses.  The  accuracy  of  the  radiation-corrected 
temperatures  is  expected  to  be  ±  25  K.  A  brief  description  of  the  experimental  procedure  is 
given  here. 

In  the  extinction  experiments  the  temperature  of  the  fuel  stream,  T\  =  378  K,  and  the 
temperature  of  the  oxidizer  stream,  X2  =  298  K.  The  distance  between  the  fuel  boundary 
and  the  oxidizer  boundary  was  L  =  10  mm.  At  some  selected  value  of  Ypp  the  flame  was 
stabilized  at  a 2  <  <32, e.  The  strain  rate  was  increased  by  increasing  V\  and  V2  until  extinction 
was  observed.  Experimental  results  are  shown  later.  Previous  autoignition  experiments  were 
conducted  with  the  mole  fraction  of  prevaporized  fuel  maintained  at  0.15  [11].  The  temperature 
at  the  fuel  boundary,  T\  =  378  K.  The  distance  between  the  fuel  boundary  and  the  oxidizer 
boundary  was  L  =  12  mm.  At  a  given  strain  rate  and  oxidizer  temperature  T2  <  X2,/  the  flow 
field  was  established.  The  temperature  at  the  oxidizer  boundary  was  gradually  increased  until 
autoignition  took  place.  Experimental  results  are  shown  later. 

Detailed  Chemical-Kinetic  Mechanism 

The  detailed  chemical-kinetic  mechanism  for  toluene  was  assembled  by  adding  the  toluene 
and  benzene  reaction  mechanism  of  Zhong  et  al.  [2-5]  to  the  C1-C4  mechanism  of  Refs.  [15] 
and  [16].  The  detailed  chemical-kinetic  mechanism  is  made  up  349  species  and  1631  reversible 
reactions.  The  entire  mechanism  will  be  available  from  our  Web  page  (see  Ref.  [17]).  The  ulti¬ 
mate  objective  of  this  work  is  to  add  toluene  as  a  fuel  component  to  detailed  chemical-kinetic 
mechanisms  for  alkanes  that  have  been  developed  at  Lawrence  Livermore  National  Labora¬ 
tories  (LLNL)  [18].  This  combined  alkane-toluene  mechanism  will  provide  a  more  complete 
mechanism  for  modeling  combustion  of  gasoline  or  diesel  fuels. 

The  toluene  mechanism  was  improved  in  several  ways.  Rate  constants  of  key  reactions  were 
estimated  using  Quantum  RRK  analysis  to  obtain  k(E)  and  master  equation  analysis  [19]  to 
evaluate  pressure  fall-off.  Reactions  analyzed  included  toluene  decomposition  reactions  and 
the  reaction  of  benzyl  radicals  with  O  atoms.  The  unimolecular  decomposition  of  toluene 
giving  benzyl  +  H  was  found  to  be  an  important  reaction  controlling  shock  tube  ignition. 
The  reaction  of  benzyl  radicals  with  O  atoms  to  give  benzyaldehyde  +  H  is  the  primary  path 
consuming  benzyl  under  shock  tube  and  flow  reactor  conditions. 

The  reaction  rate  constants  for  phenyl  +  O2  system  developed  by  Bozzelli  et  al.  [19]  were 
incorporated.  These  reaction  channels  included  new  products  for  which  consumption  reac¬ 
tions  were  added.  The  phenyl  +  O2  rate  constants  were  found  to  play  an  important  role 
under  counterflow-ignition  conditions. 
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The  rate  constant  for  the  initiation  reaction  (toluene  +  O2  #  benzyl  +  HO2)  was  updated 
(k  =  9.3xl08  T1-301  exp(— 40939.0 cal/RT)  cm3/(mol-s))  using  one  half  the  rate  constant  of 
allylic  isobutene  +  O2  [20].  Our  estimate  is  50  percent  higher  than  the  estimate  of  Walker  [21] 
at  1000  K.  Reactions  to  consume  bicyclopentadiene  were  also  added. 

Comparisons  between  Numerical  Calculations  and  Experiments 

Shock  Tube  Comparisons 

The  ignition  in  the  shock  tube  was  simulated  using  the  Senkin  code  [22]  assuming  constant 
volume  combustion  after  the  reflected  shock  passes  through  the  mixture.  The  measured  results 
are  compared  to  predicted  results  in  Fig.  2.  In  the  experiments,  10  percent  of  the  maximum 
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Figure  2:  Predicted  (lines  and  filled  symbols)  and  measured  (unfilled  symbols)  ignition  delay 
times  of  toluene/02/Ar  mixtures  under  shock  tube  conditions.  (T  =  reflected  shock  temper¬ 
ature  in  Kelvin). 

OH  emission  was  used  as  an  indication  of  ignition.  The  computed  ignition  delay  time  was 
obtained  using  the  inflection  in  the  temperature  profile  as  an  indication  of  ignition.  When 
10  percent  of  the  OH  maximum  concentration  was  used  in  the  calculations,  almost  identical 
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results  were  obtained  over  the  temperature  range  considered.  We  did  not  attempt  to  simulate 
OH  emission.  The  predicted  ignition  delay  times  compared  reasonably  well  with  experimental 
values  for  the  fuel-rich  case,  but  were  too  long  for  the  fuel-lean  case. 

Standard  Senkin  sensitivity  analysis  [22]  was  performed  to  determine  the  reaction  rate 
constants  that  control  the  oxidation  of  toluene  under  shock  tube  conditions.  The  Senkin  code 
gives  the  change  in  species  concentration  for  an  incremental  change  in  reaction  rate  constant. 
The  OH  radical  concentration  was  chosen  as  an  indicator  of  the  overall  reactivity  of  the  system. 
An  alternate  choice  of  toluene  concentration  produced  similar  results.  The  sensitivity  of  OH 
concentration  to  a  change  in  reaction  rate  constant  is  shown  in  in  Fig.  3.  The  results  are 


H  +  02  =  OH  +  O 
C6H5CH3  +  H  =  C6H5CH2  +  H2 
C6H5CH3  =  C6H5CH2  +  H 
C6H5CHO  =  C6H5CO  +  H 
C6H5CH3  +  OH  =  C6H5CH2  +  H20 
C6H5CH3  +  02  =  C6H5CH2  +  H02 
C6H5CH2  +  H02  =  C6H5CH20  +  OH 
C6H5CH2  +  02  =  C6H5CHO  +  OH 
C6H5CH3  =  C6H5  +  CH3 
C6H5CHO  +  OH  =  C6H5CO  +  H20 
C6H5CH2.  =  CY-C5H5  +  C2H2 
HCO  +  02  =  CO  +  H02 
C6H5CH3  +  H  =  C6H6  +  CH3 
HCO  +  M  =  H  +  CO  +  M 
CY-C5H5+02<=>C=CC=C=0+HC0 
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Figure  3:  Sensitivity  of  the  OH  radical  concentration  to  changes  in  individual  rate  constants 
under  the  shock  tube  conditions  shown  in  Fig.  2  with  1000/T  =  0.67,  where  T  is  given  in 
Kelvin. 

shown  when  the  toluene  is  about  half  consumed.  Positive  sensitivities  indicate,  an  increase  in 
rate  constant  increases  the  OH  concentration  and  accelerates  the  overall  rate  of  reaction,  while 
negative  sensitivities  indicate  an  increase  in  rate  constant  decreases  the  OH  concentration  and 
retards  the  overall  rate  of  reaction.  The  reaction  giving  the  highest  sensitivity  is  H  +  O2  chain 
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branching  (Fig.  3).  The  sensitivities  of  the  other  controlling  reactions  can  be  understood  in 
relationship  to  it.  The  abstraction  reaction 


c6h5ch3  +  h  =  c6h5ch2-  +  h2, 

removes  H  atoms  that  would  otherwise  lead  to  chain  branching  via  H  +  02  so  it  exhibits 
a  high  negative  sensitivity  and  retards  ignition.  We  used  the  rate  constant  from  Baulch  et 
al.  [23].  This  rate  constant  has  an  uncertainty  of  about  a  factor  of  2  at  1700  K.  The  toluene 
decomposition  channel, 


C6H5CH3  =  C6H5CH2-  +  H, 

produces  H  atoms  that  can  then  lead  to  chain  branching  and  thus  exhibits  a  high  positive 
sensitivity  and  accelerates  ignition.  This  reaction  was  analyzed  carefully  with  quantum  RRK 
and  master  equation  treatment  because  of  its  importance  under  shock  tube  conditions.  The 
reaction  C6H5CHO  =  CgHsCO-  +  H  also  produces  H  atoms  and  therefore  exhibits  a  signifi¬ 
cant  positive  sensitivity.  The  inclusion  of  this  reaction  significantly  improved  the  shock  tube 
predictions  at  high  temperature.  The  reaction  of  toluene  with  OH  shows  a  large  negative 
sensitivity  and  thus  retards  ignition.  At  fuel-lean  conditions,  the  reaction  CgH5CH3  +  02 
=  C6H5CH2-  +  H02  shows  a  high  positive  sensitivity.  Further  theoretical  or  experimental 
investigations  of  this  reaction  and  the  reaction  of  toluene  with  H  atoms  may  help  improve 
agreement  with  the  shock  tube  experiments  under  fuel-lean  conditions. 

Flow  Reactor  Comparisons 

The  flow  reactor  was  simulated  using  the  Senkin  code  [22] .  The  predicted  results  are  compared 
to  the  measurements  of  Klotz  et  al.  [1]  in  Figs.  4  and  5.  The  results  of  the  fuel  profile  look 
reasonable  and  the  early  appearance  of  the  benzylaldehyde  peak  compared  to  the  benzene 
peak  is  predicted.  Many  of  the  peak  intermediate  concentrations  are  predicted  within  a  factor 
of  two  of  the  measurements,  except  for  benzyladehyde  which  is  within  a  factor  of  three  and 
1,3-butadiene  which  is  predicted  to  be  in  very  low  concentrations  compared  to  the  measure¬ 
ments. 

Senkin  sensitivity  analysis  [22]  was  performed  to  determine  the  reaction  rate  constants  that 
control  the  oxidation  of  toluene  under  flow  reactor  conditions.  The  toluene  concentration  was 
chosen  as  an  indicator  of  the  overall  reactivity  of  the  system.  The  sensitivity  of  the  toluene 
concentration  at  135  ms  to  a  change  in  reaction  rate  constant  is  shown  in  Fig.  6.  At  135  ms,  the 
initial  toluene  is  about  48  percent  consumed.  Negative  sensitivities  indicate  that  the  reaction 
accelerates  the  overall  rate  of  reaction  and  positive  sensitivities  indicate  the  opposite  effect. 
The  reaction  exhibiting  the  highest  sensitivity  is  again  the  H  +  02  chain  branching  reaction. 
The  second  most  sensitive  reaction  is  the  reaction  of  benzyl  and  02: 
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Figure  4:  Predicted  and  measured  concentration  histories  under  flow  reactor  conditions.  (</>  = 
0.76,  N2  98%,  initial  temperature  =  1173  K,  atmospheric  pressure,  time  is  residence  time  in 
flow  reactor).  The  symbols  represent  measurements  of  Klotz  et  al.  [1].  The  lines  are  results 
of  numerical  calculations. 
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Figure  5:  Predicted  and  measured  concentration  histories  under  flow  reactor  conditions.  (</>  = 
0.76,  N2  98%,  initial  temperature  =  1173  K,  atmospheric  pressure,  time  is  residence  time  in 
flow  reactor).  The  symbols  represent  measurements  of  Klotz  et  al.  [1].  The  lines  are  results 
of  numerical  calculations. 
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H+02=0H+0 

C6H5CH2.+02=C6H50.+CH20 

HC0+02=C0+H02 

HCO+M=H+CO+M 

C6H5CH3+H=C6H5CH2+H2 

C6H5CH0+0H=C6H5C0+H20 

C6H5CH3+02=C6H5CH2+H02 

C6H5CH3+0H=C6H5CH2+H20 

C6H5CH3+H=C6H6+CH3 

C6H5CH2.+0=C6H5CH0+H 

CH3+02=CH20+0H 

C6H5CH3=C6H5CH2+H 

C6H5CH2+H02=C6H5CH20+0H 

C6H5CHO=C6H5CO+H 

CH3+CH3(+M)=C2H6(+M) 

C6H5+02=C6H50.+0 

C6H5CHO+H=H2+C6H5C.O 

C6H5+02<=>0=C.C=CC=CC=0 

C6H5CH2.+0H=C6H5CH20H 
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Figure  6:  Sensitivity  of  the  toluene  concentration  to  changes  in  individual  rate  constants 
under  flow  reactor  conditions,  (equivalence  ratio  =  0.76,  N2  98%,  initial  temperature  = 
1173K,  atmospheric  pressure,  residence  time  is  135  ms). 
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c6h5ch2-  +  o2  =  c6h5o-  +  ch2o. 

The  rate  constant  of  this  reaction,  k  =  5.30x  1013T_107exp(— 1 0840 cal /-RT),  was  taken  from 
Zhong  [4],  This  reaction  occurs  through  a  4-membered  transition  state  where  the  terminal 
oxygen  on  the  benzyl  peroxy  radical  adds  to  the  benzene  ring. 

The  HCO  +  02  =  CO  +  H02  reaction  promotes  fuel  oxidation  because  it  produces  H02 
radicals  that  react  with  benzyl  radicals  via, 

C6H5CH2-  +  H02  =  C6H5CH20-  +  OH 

a  reaction  with  a  negative  sensitivity  that  also  promotes  fuel  oxidation.  The  reaction  of  toluene 
with  H  atoms  exhibits  a  positive  sensitivity  and  retards  the  overall  fuel  oxidation  rate,  as  seen 
under  shock  tube  conditions. 

Comparisons  with  Strained  Laminar  Flow  under  Nonpremixed  Conditions 

The  detailed  chemical-kinetic  mechanism  was  tested  by  comparing  computed  results  with  ex¬ 
periments  performed  under  strained,  nonpremixed  conditions.  Numerical  calculations  were 
carried  out  using  FlameMaster  [24].  The  formulation  of  the  numerical  problem  is  summa¬ 
rized  elsewhere  [24,25].  At  both  ends  of  the  computational  domain,  the  mass  fractions  of 
the  reactants  and  the  normal  components  of  the  flow  velocity  were  specified.  The  values  of 
the  tangential  component  of  the  flow  velocity  at  both  ends  were  set  equal  to  zero  (the  so- 
called  plug- flow  boundary  conditions).  The  characteristic  strain  rate  at  the  stagnation  plane 
was  calculated  using  Eq.  (1).  Solutions  could  not  be  obtained  with  the  349  species  detailed 
chemical-kinetic  mechanism  due  to  numerical  ” stiffness”  problems.  The  detailed  chemical  ki¬ 
netic  mechanism  was  simplified  using  the  NIST  XSenkplot  [26].  Simplified  mechanisms  were 
obtained  under  shock  tube  conditions  and  flow  reactor  conditions  and  combined  to  yield  a  58 
species  mechanism  made  up  of  221  reversible  reactions.  Any  reaction  in  the  detailed  mecha¬ 
nism  that  involved  only  these  58  species  was  retained  in  the  simplified  mechanism.  Ignition 
delay  times  were  calculated  using  the  simplified  mechanism  and  the  detailed  mechanism  at 
conditions  similar  to  those  employed  in  the  shock-tube  experiments.  At  all  conditions  the 
differences  were  found  to  be  less  than  14%. 

Figure  7  shows  the  mass  fraction  of  toluene  in  the  fuel  stream  at  extinction,  Yp  i ,  as  a  func¬ 
tion  of  the  strain  rate,  a2;e.  The  symbols  represent  measurements.  The  solid  line  represents 
results  of  numerical  calculations  performed  using  the  simplified  chemical-kinetic  mechanism. 
At  a  given  value  of  Yfa  the  calculated  value  of  a2je  is  lower  than  that  measured.  Figure  8 
shows  the  oxidizer  temperature  at  autoignition,  as  a  function  of  the  strain  rate,  a2y.  The 
symbols  represent  measurements  reproduced  from  Refs.  [11,  27].  The  solid  line  represents  re¬ 
sults  of  numerical  calculations  performed  using  the  simplified  chemical-kinetic  mechanism.  At 
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Figure  7:  Mass  fraction  of  toluene  in  the  fuel  stream  at  extinction,  1f,i>  as  a  function  of  the 
strain  rate,  <32, e-  The  symbols  represent  measurements.  The  solid  line  represents  results  of 
numerical  calculations  performed  using  the  simplified  chemical-kinetic  mechanism  made  up  of 
58  species. 


Strain  Rate  at  Ignition,  a1 ,  [s  '] 


Figure  8:  Oxidizer  temperature  at  autoignition,  T2J,  as  a  function  of  the  strain  rate,  02,/. 
The  symbols  represent  measurements  [11,27].  The  solid  line  represents  results  of  numerical 
calculations  performed  using  the  simplified  chemical-kinetic  mechanism  made  up  of  58  species. 
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a  given  value  of  the  oxidizer  temperature  the  calculated  value  of  a2j  is  lower  than  that  mea¬ 
sured.  Thus  in  both  extinction  and  ignition  experiments,  the  numerical  model  predicts  that 
toluene/air  is  overall  less  reactive  than  observed  in  the  experiments.  This  result  is  consistent 
with  comparisons  between  the  model  predictions  and  experimental  results  from  the  shock  tube 
above  1400K. 

The  sensitivity  results  for  autoignition  under  nonpremixed  conditions  are  given  in  Fig.  9. 
The  analysis  performed  by  FlameMaster  [24]  gives  the  change  in  maximum  OH  concentration 


H  +  02  =>  O  +  OH 
C6H5CH2.+H02  =>  C6H5CH20.+0H 
H02  +  0H  =>  H20  +  02  j 
C6H5CHO+  OH  =>  H20+C6H5C0: 

H02  +  O  =>  OH  +  02  - 
HCO  +  M  =>  H  +  CO  +  M 
C6H5  +  02  =>  C6H50.  +  O  1 
HCO  +  02  =>  CO  +  H02 
C6H5+02  =>  0=C.C=CC=CC=0: 
C6H5CH2.  +  OH  =  C6H5CH20H  - 
H  +  02  +  M  =>  H02  +  M 
C6H5CH3  +  H  =>  C6H5CH2.  +  H2 
C6H5CH3  +  O  =>  C6H5CH2.  +  OH 
C6H5CH3  +  OH  =  C6H5CH2.  +  H20: 
C6H5CHO  =>  C6H5C.O  +  H: 
C6H5CH3  +  H  =>  C6H6  +  CH3  - 


-10-5  0 


10  15 


Sensitivity 


Figure  9:  Sensitivity  of  the  OH  radical  concentration  to  changes  in  individual  rate  constants 
under  nonpremixed  conditions  near  autoignition  (strain  rate  =  400  s-1,  oxidizer  temperature 
T2  =  1357  K). 

for  an  incremental  change  in  rate  constant  as  indicated  in  Fig.  9.  In  FlameMaster,  forward  and 
reverse  rate  constants  are  considered  as  separate  parameters  for  the  purposes  of  sensitivity 
analysis.  The  analysis  was  performed  for  a  reactive  flow  solution  very  near  autoignition. 
Positive  sensitivities  indicate  an  increase  in  rate  constant  increases  the  OH  concentration  and 
accelerates  the  overall  rate  of  reaction,  while  negative  sensitivities  indicate  an  increase  in  rate 
constant  decreases  the  OH  concentration  and  retards  the  overall  rate  of  reaction.  The  H  + 
O2  chain  branching  reaction  again  gives  very  high  sensitivity  as  under  shock  tube  and  flow 
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reactor  conditions.  The  reaction  benzyl  radicals  with  HO2, 

C6H5CH2-  +  H02  =  C6H5CH20-  +  OH 

exhibits  a  particularly  high  sensitivity  under  counterflow  ignition  conditions.  The  reaction  of 
benzaldehyde  with  OH  shows  a  high  positive  sensitivity  and  promotes  ignition. 

The  reactions  that  give  significant  sensitivity  under  nonpremixed  conditions  but  not  under 
flow  reactor  or  shock  tube  conditions  are  reactions  that  involve  the  phenyl  +  02  reaction: 

c6h5  +  o2  ^  c6h5o-  +  o, 
c6h5  +  o2  ^  o=c-c=cc=cc=o. 

An  important  uncertainty  in  the  phenyl  +  02  system  is  the  branching  ratio  between  the  two 
product  channels  given  above.  The  0-atom  product  channel  is  accelerating  and  the  other 
ring-opening  product  channel  is  retarding  (Fig.  9).  The  branching  ratio  depends  on  the 
relative  barrier  heights  for  the  phenylperoxy  radical  going  to  0-0  bond  breakage  versus  to 
ring  opening  [19]. 

Conclusions 

The  detailed  chemical-kinetic  model  for  toluene  was  improved  by  adding  a  more  accurate  de¬ 
scription  of  the  phenyl  +  02  reactions,  toluene  decomposition  reactions  and  the  benzyl  radical 
+  O  reaction.  These  reactions  have  been  analyzed  by  quantum  RRK  to  obtain  k(E)  and  by 
master  equation  for  fall-off.  New  data  obtained  under  nonpremixed  conditions  in  a  counterflow 
configuration  and  obtained  in  a  shock  tube  were  compared  to  results  of  the  detailed  chemical- 
kinetic  model.  Sensitivity  analysis  was  used  to  identify  reactions  whose  rate  constants  control 
the  overall  rate  of  oxidation.  This  information  can  be  used  to  obtain  future  improvements  of 
the  agreement  between  the  model  and  experiments.  The  reactions  that  exhibited  high  sen¬ 
sitivity  (besides  H  +  02  chain  branching)  included  toluene  decomposition  reactions,  toluene 
reaction  with  H,  OH  and  02,  and  benzyl  reaction  with  H02  and  02.  Also,  phenyl  +  02 
reactions  were  shown  to  be  important  during  ignition  in  a  nonpremixed,  counterflow  system. 
Although  the  benzyl  +  O  reaction  did  not  appear  in  the  sensitivity  analysis,  it  was  primary 
reaction  consuming  benzyl  under  shock  tube  and  flow  reactor  conditions. 
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Appendix  B:  The  Structure  of  Nonpremixed  n- Decane  Flames 

K.  Tanoue,  R.  Seiser,  K.  Seshadri 

Abstract 

Experiments  are  conducted  on  flames  stabilized  between  two  countcrflowing  streams. 

The  fuel  stream  is  a  mixture  of  prevaporized  n-decane  and  nitrogen,  and  the  oxidizer 
stream  is  a  mixture  of  air  and  nitrogen.  Concentration  profiles  of  C10H22,  O2,  N2,  CO2, 

CO,  H20,  CO,  CH4,  and  hydrocarbons  ranging  from  C2  up  to  Cg,  are  measured.  The  mea¬ 
surements  are  made  by  removing  gas  samples  from  the  flame  using  a  quartz  microprobe 
and  analyzing  the  samples  using  a  gas  chromatograph.  Temperature  profiles  are  mea¬ 
sured  using  a  thermocouple.  Critical  conditions  of  extinction  and  critical  conditions  of 
autoignition  are  measured.  Numerical  calculations  are  performed  using  detailed  chemistry 
to  determine  the  flame  structure  and  to  obtain  values  for  extinction  and  autoignition. 

Introduction 

The  combustion  of  n-decane  has  been  studied  experimentally  by  various  investigators.  Previ¬ 
ous  studies  were  conducted  on  shock-tubes  [1],  jet-stirred-reactors  [2],  flow-reactors  [3],  burner 
stabilized  [4,  5],  and  freely  propagating  premixed  flames  [6].  The  present  experimental  investi¬ 
gation  reports  experimental  data  obtained  for  a  nonpremixed  n-decane  flame  in  a  counterflow 
configuration.  In  the  following,  the  experimental  method  is  described,  and  data  for  stable 
species  and  temperature  are  reported.  In  addition,  data  for  extinction  and  autoignition  are 
provided. 

Description  of  Experimental  and  Numerical  Studies 

Studies  are  conducted  on  nonpremixed  n-decane  flames  stabilized  in  the  counterflow  config¬ 
uration.  For  temperature  and  species  profiles  the  fuel  stream  contains  n-decane  with  a  mole 
fraction  of  Xc10h22= 0-08  and  nitrogen.  The  oxidizer  stream  is  air.  The  strain  rate  02=100 s-1 
and  is  calculated  by  the  formula  [7] 


a  2 


21W  \Vi\Vp!\ 
L  \  |E2|  yW 


(1) 


Here,  V  denotes  the  velocity  and  p  the  density.  Subscripts  1  and  2  refer  to  the  fuel  stream  and 
oxidizer  stream,  respectively.  The  separation  distance  between  the  two  ducts  is  L=10mm. 
The  momenta  of  the  two  opposing  streams  are  balanced  as  described  by  p\  Vf  =  P2V2  ■  Plug- 
flow  conditions  are  assumed.  Details  and  dimensions  of  the  counterflow  burner  are  given  in 
Refs.  [8]  and  [9].  Figure  1  shows  the  setup  for  temperature  and  species  profile  measurements. 
Concentrations  of  stable  species  are  measured  by  removing  gas  samples  from  the  reaction 
zone  using  a  quartz  microprobe  and  analyzing  them  in  a  gas  chromatograph.  The  microprobe 
has  a  tip  with  an  outer  diameter  of  166  microns  and  an  inner  diameter  of  83  microns.  The 
tip  is  placed  at  a  location  of  5  mm  off  the  centerline  of  the  ducts  to  avoid  disturbance  of 
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Figure  1:  Temperature  and  species  profile  measurements  in  a  nonpremixed  n-decane  flame. 

the  flow-field  and  influencing  the  reactive  diffusive  balance  in  the  vicinity  of  the  probe.  The 
species  concentrations  and  temperatures  are  assumed  to  be  only  functions  of  the  coordinate 
along  the  centerline  [10].  The  samples  drawn  from  the  flame  are  quantified  using  a  SRI 
8610C  gas  chromatograph.  The  instrument  is  equipped  with  with  a  4.5ft  mole-sieve  (80/100 
mesh)  -  for  separating  H2,  O2  and  Ar,  N2,  CH4,  and  CO  -  and  a  12ft  Porapaq  Q  column  for 
separating  all  other  species.  Temperature  programming  and  valve  switching  is  employed  to 
optimize  the  separation  performance  of  both  columns.  The  species  eluding  from  the  column  are 
quantified  using  a  thermo-conductivity-detector  (TCD)  and  a  flame-ionization-detector  (FID). 
The  chromatograms  are  analyzed  using  self  written  software,  and  by  comparing  with  runs  of 
known  samples,  the  absolute  mole  fractions  are  determined.  For  species  where  no  calibration 
gases  were  used,  we  determined  the  calibration  factors  by  inter-  or  extra-polating  around 
known  species.  This  was  done  by  plotting  the  FID-calibration  factors  as  a  function  of  the 
number  of  C-atoms.  In  one  curve  for  the  paraffins,  species  measured  comprise  ethane,  propane, 
n-heptane,  n-octane,  and  n-decane.  A  second  curve  for  the  alkenes  was  similarly  established 
by  calibrating  ethene,  propene,  1-butene,  and  1-hexene.  All  species  of  C4-hydrocarbons  are 
reported  as  one  concentration,  because  individual  separation  and  calibration  was  not  possible. 
Similarly  C5-,  C6-,  C7-,  Cs-,  and  Cg-hydrocarbons  are  reported,  here  assuming  that  they 
are  mostly  comprised  by  1-alkenes,  and  thus  using  the  1-alkene  calibrations.  Ethene  and 
ethyne  appeared  as  single  peak  and  are  reported  as  their  sum  using  the  calibration  for  ethene. 
Previous  analysis  using  a  mass-spectrometer  has  shown  that  for  a  n-heptane  flame  the  amount 
of  ethene  is  slightly  larger  than  that  of  ethyne  [11].  Argon  eluded  together  with  oxygen  from 
the  columns.  By  assuming  the  Lewis  number  of  argon  and  the  Lewis  numbers  of  the  main 
species  to  be  equal,  the  amount  of  argon  at  any  point  in  the  flow-field  was  determined  from 
the  mixture-fraction  (calculated  using  the  carbon  mass  fraction  from  the  measured  species) 
and  the  amount  of  argon  in  the  oxidizer  stream  (Yat =0.01313).  Water  appeared  as  peak 
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with  a  long  tail  on  the  TCD-detector.  The  area  of  this  peak  was  determined  by  using  a 
linear  baseline  starting  at  the  onset  of  the  water  signal  and  parallel  to  the  baseline  for  a 
measurement  containing  no  water.  By  choosing  a  proper  balance  of  pressure  and  temperature 
of  the  sampling  probe,  sampling  line,  and  sample  loop,  we  were  able  to  measure  all  water 
without  experiencing  condensation.  This,  however,  was  not  possible  for  certain  measurements 
containing  n-decane,  as  its  boiling  point  is  much  higher  than  that  of  water.  We  therefore  only 
report  values  of  n-decane  smaller  than  Xc10h22= 0-03.  Closer  to  the  fuel-side,  the  amounts 
of  n-decane  are  estimated  using  a  numerical  calculation  of  a  Burke-Schumann  (subscript  BS) 
solution.  In  these  locations  the  amounts  of  all  other  species  can  then  be  determined  by  using 
the  measured  values  and  normalizing  them  so  that  their  sum  equals  (1.0  -  Xc10h22,bs)-  The 
expected  accuracy  for  the  maximum  concentrations  of  most  species  is  expected  to  be  better 
than  ±10%.  This  includes  species  that  can  be  clearly  identified  on  either  the  FID  or  TCD- 
detector.  Hydrogen  gives  a  very  small  signal  on  the  TCD-detector  and  it  is  estimated  to  be 
accurate  to  ±25%.  The  expected  accuracy  for  water  is  asurned  to  be  ±20%.  At  this  point 
the  accuracy  for  C7-C9  species  is  not  quantifiable,  as  their  absolute  concentrations  are  small 
and  possibly  influenced  by  radical  recombination  in  the  probe.  The  temperature  profile  was 
measured  using  an  uncoated  Pt-Ptl3%Rh  (Type  R)  thermocouple.  The  wire  diameter  is 
25  microns  and  the  bead  diameter  is  90  microns.  The  measurement  is  corrected  for  radiation 
by  assuming  spherical  shape  of  the  bead  and  a  constant  emissivity  of  0.2.  Catalytic  effects  are 
neglected.  The  absolute  accuracy  of  the  temperature  measurement  is  expected  to  be  better 
than  ±80  K.  The  location  of  the  sampling  probe  and  the  thermocouple  in  the  flow-field  is 
determined  using  a  digital  photo  camera.  The  size  of  one  pixel  corresponds  to  a  distance  in 
the  flow- field  of  approximately  17  microns. 

Extinction  measurements  are  performed  for  various  n-decane  concentrations  in  the  fuel 
stream.  The  oxidizer  stream  is  air.  For  a  certain  composition  of  the  reactant  streams,  a 
flame  is  established,  and  the  the  strain  rate  is  increased  until  extinction  is  observed.  The 
autoignition  experiments  are  performed  for  a  mole  fraction  of  n-decane  Xc10#22=0.15.  The 
oxidizer  stream  is  air,  and  it  is  heated  until  autoignition  is  observed.  In  the  extinction  and 
ignition  experiments  the  separation  distance  of  the  ducts  is  10  mm  and  12  mm,  respectively. 
The  temperature  of  the  fuel  stream  is  400  K  in  the  extinction  experiments,  and  408  K  in  the 
ignition  experiments. 

Numerical  calculations  are  performed  using  detailed  chemistry,  applying  the  same  bound¬ 
ary  conditions  as  in  the  experiments.  For  the  computations  the  computer  program  FlameMas- 
ter  is  used,  which  was  developed  at  RWTH- Aachen  [12].  The  conservation  equations  of  mass, 
momentum  and  energy  and  the  species  balance  equations  used  in  the  formulation  of  the  nu¬ 
merical  problem  are  summarized  elsewhere  [12-14],  The  species  balance  equations  include 
thermal  diffusion  and  the  energy  conservation  equation  includes  radiative  heat  losses  from 
carbon  dioxide  and  water  vapor  [12].  Plug-flow  boundary  conditions  are  used.  Buoyancy  is 
neglected.  As  chemical-kinetic  mechanism,  a  detailed  mechanism  by  Bikas  et  al.  [15]  contain¬ 
ing  347  reversible  reactions  among  66  species,  is  used. 


35 


Results  and  Discussion 

Figure  2  shows  experimental  data  on  the  temperature  field  in  the  nonpremixed  n-decane  flame. 


Figure  2:  Temperature  measurement  in  a  nonpremixed  n-decane  flame.  The  triangles  represent 
the  raw-data  of  the  measured  temperature,  the  spheres  represent  the  temperature  corrected 
for  radiation.  The  lines  are  results  of  numerical  calculations  using  detailed  chemistry. 

The  triangles  show  the  temperature  of  the  thermocouple  as  read  from  the  temperature  display 
unit.  After  a  radiation  correction  is  applied,  that  takes  into  account  the  radiation  from  the 
surface  of  the  thermocouple,  the  gas  temperature  is  obtained,  shown  as  spheres  in  Fig.  2. 
The  corrected  temperature  agrees  well  with  the  computed  one,  concerning  the  maximum 
temperature,  shift  and  broadness  of  the  temperature  profile.  Figure  3  shows  experimental  data 
on  the  main  species  in  the  n-decane  flame.  The  agreement  between  measured  and  calculated 
concentrations  is  very  good.  CO2,  H2O,  and  CO  are  the  main  products  of  combustion  and  are 
usually  well  predicted  by  chemical-kinetic  mechanisms.  The  comparison  in  Fig.  3  is  therefore 
an  indication  of  the  flow-field  and  the  influence  by  the  microprobe.  The  influence  of  the 
microprobe  on  the  reactive-diffusive  balance  is  usually  high  closer  to  the  ducts  where  the  axial 
velocity  -  perpendicular  to  the  probe  -  is  higher.  In  the  reaction  zone,  heat  loss  to  the  probe 
can  cause  local  quenching  and  cause  stronger  oxygen  leakage.  From  the  profile  for  oxygen 
it  can  be  seen  that  the  measured  leakage  is  slightly  higher  than  predicted  by  the  numerical 
calculation.  The  chemical-kinetic  mechanism,  however,  predicts  a  higher  extinction  strain 
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Figure  3:  Experimental  data  showing  the  mole  fraction  of  C10H22,  O2,  CO2,  H2O,  and  CO 
as  a  function  of  distance.  The  points  represent  measurements  and  the  lines  are  results  of 
numerical  calculations  using  the  chemical-kinetic  mechanism  by  Bikas  et  al.  [15] . 

rate,  and  therefore  less  leakage  at  the  measured  strain  rate  of  100  s_1.  The  influence  of  the 
microprobe  can  therefore  be  assumed  to  be  minimal.  Figures  4  and  5  show  experimental 
data  on  intermediate  species  of  small  molecular  weight.  The  profiles  of  C2H2+C2H4,  H2, 
and  C3H8  agree  well  with  numerical  calculations.  The  chemical-kinetic  mechanism  under¬ 
predicts  CH4,  C2H6,  and  CsHq  by  about  30-70%.  Figures  6  and  7  show  experimental  data  on 
concentrations  of  larger  stable  hydrocarbon  species.  Radical  species  are  usually  small  and  will 
recombine  in  the  sampling  probe  to  form  stable  species.  The  numerical  computations  for  C4- 
species  somewhat  over-predict  measured  concentrations,  while  computations  of  stable  C5-  and 
Cg-species  are  lower  than  the  experimental  data.  Measured  stable  C7-hydrocarbons  are  much 
larger  than  predicted.  The  calculated  values  for  the  sum  of  all  C7-species  (including  radical 
species)  are  everywhere  below  a  mole  fraction  of  10-8.  For  Cs-  and  Cg-species,  no  stable 
species  appeared  in  the  chemical-kinetic  mechanism.  Figure  8  shows  experimental  data  on 
extinction  of  a  nonpremixed  n-decane  flame  with  variable  fuel  concentrations.  The  results  are 
compared  with  nunrercal  calculations  using  the  chemical-kinetic  mechanism  of  Bikas  et  al.  [15]. 
The  mechanism  agrees  reasonably  well  with  experiments.  Figure  9  shows  experimental  data 
on  autoignition  of  a  nonpremixed  n-decane  /  air  system.  The  numerical  calculations  agree 
well  with  experiments,  especially  at  lower  strain  rates.  At  higher  strain  rates  the  mechanism 
predicts  slightly  lower  autoignition  temperatures. 


37 


Figure  4:  Experimental  data  showing  the  mole  fraction  of  C2H2+C2H4,  CH4,  and  C2H6  as  a 
function  of  distance.  The  points  represent  measurements  and  the  lines  are  results  of  numerical 
calculations. 


Figure  5:  Experimental  data  showing  the  mole  fraction  of  H2,  CsHg,  and  C3H8  as  a  func¬ 
tion  of  distance.  The  points  represent  measurements  and  the  lines  are  results  of  numerical 
calculations. 
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Figure  6:  Experimental  data  showing  the  mole  fraction  of  C4-,  C5-,  and  C6-hydrocarbons  as  a 
function  of  distance.  The  points  represent  measurements  and  the  lines  are  results  of  numerical 
calculations.  The  numerically  computed  temperature  is  shown  for  reference. 
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Figure  7:  Experimental  data  showing  the  mole  fraction  of  C7-,  Cs-,  and  Cg-hydrocarbons 
as  a  function  of  distance.  The  points  represent  measurements.  The  numerically  computed 
temperature  is  shown  for  reference. 
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Figure  8:  Extinction  of  nonpremixed  n-decane  flame.  The  mole  fraction  of  n-decane  in  the  fuel 
stream  is  plotted  versus  the  strain  rate.  The  oxidizer  stream  is  air.  The  temperatures  of  the 
fuel  stream  and  oxidizer  stream  are  400  K  and  298  K.  The  symbols  represent  measurements 
and  the  lines  are  results  of  numerical  calculations. 
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Appendix  C:  Autoignition  and  Extinction  of  Hydrocarbon  Fuels  in  Nonpremixed  Systems 

S.  Humer,  R.  Seiser,  K.  Seshadri 

Abstract 

Experimental  studies  are  carried  out  on  extinction  and  autoignition  of  liquid  hydrocarbon  fuels  in  nonpremixed  systems. 
The  counterflow  configuration  is  employed  in  the  experimental  studies.  In  this  configuration  an  oxidizer  made  up  of  air 
and  nitrogen  flows  over  the  surface  of  a  pool  of  liquid  fuel.  Critical  conditions  of  extinction  and  autoignition  are 
measured  for  n-heptane,  iso-octane,  «-octane,  « -decane,  n-dodecane,  w-hexadecane,  o -xylene,  JP-10,  diesel,  and  JP-8. 

Introduction 

Fundamental  studies  on  extinction  and  autoignition  of  strained  flames  provide  knowledge  for  modeling  turbulent 
combustion.  Previous  studies  on  autoignition  of  reactant  mixtures  were  focused  on  measuring  ignition  delay  time  in 
shock  tubes  [1]  and  rapid  compression  machines  P].  These  studies  do  not  characterize  the  influence  of  strain  on 
autoignition.  In  this  study  the  influence  of  strain  on  the  critical  conditions  of  extinction  and  autoignition  of  flames  is 
investigated. 

Increasingly  diesel  is  the  fuel  of  choice  in  the  automobile  industry  in  Europe  in  view  of  fuel  economy  and  higher 
efficiency  of  diesel  engines  compared  to  gasoline  engines.  Up  to  70  %  of  the  new  registered  cars  are  equipped  with 
diesel  engines.  The  high  complexity  of  higher  hydrocarbon  fuels  makes  it  useful  to  develop  a  model  fuel  of  chemical 
characteristics  similar  to  that  of  diesel  to  keep  calculation  efforts  low.  These  model  fuels  can  be  used  in  test  engines  as 
well  as  in  numerical  simulations  of  engine  performance  and  rates  of  pollutant  formation.  The  study  should  provide 
information  about  fuels  to  be  used  as  a  base  for  a  diesel  model  fuel. 

The  researched  fuels  were  investigated  in  a  simple  geometry,  a  laminar  flame  in  a  counterflow  configuration.  A 
previous  research  was  conducted  in  a  similar  configuration  by  R.  Seiser  et  al.  [3].  Although  diesel  engines  burn  the  fuel 
in  turbulent  spray  combustion  at  high  pressure,  the  flame  on  the  surface  of  a  tiny  droplet  can  be  considered  as  a  laminar 
non-premixed  flame  [4]. 

Description  of  Experimental  Setup 

The  counterflow  configuration  is  employed  in  the  experimental  studies.  Figure  1  shows  a  schematic  illustration  of 
the  counterflow  configuration.  In  this  design  a  duct  is  mounted  above  a  liquid  pool.  The  oxidizer  stream  is  injected  from 
this  upper  duct  called  the  oxidizer  duct.  The  oxidizer  stream,  made  up  of  air  and  nitrogen,  flows  over  the  surface  of  a 
pool  of  liquid  fuel  at  atmospheric  pressure.  Concentric  ducts  surround  the  oxidizer  duct  and  the  cup  of  the  liquid  pool. 
This  outer  duct  provides  a  curtain  flow  of  nitrogen,  N2,  surrounding  the  inner  flow.  This  curtain  stream  of  N2  shields  the 
inner  flow,  containing  the  reactants,  from  influences  of  the  surrounding.  A  stagnation  plane  is  formed  and  the  exhaust 
products  are  drawn  into  the  exhaust  system  (not  shown). 


Figure  1  Schematic  illustration  of  the  liquid  pool  configuration. 
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In  the  exhaust  system  the  exhaust  gases  get  cooled  down  by  water  sprays  to  prevent  the  unburned  reactants  from  ignition. 
To  achieve  a  uniform  velocity  profile  at  the  end  of  the  oxidizer  duct,  three  layers  of  a  fine  screen  of  200  mesh  are 
mounted  at  the  duct  exit. 

The  distance  between  the  surface  of  the  liquid  pool  and  the  exit  of  the  oxidizer  duct  is  L.  At  the  exit  of  the  oxidizer 
duct,  the  value  of  the  injected  velocity  is  V2,  the  temperature  T2,  the  density  P2,  and  the  mass  fraction  of  oxygen  Yo2,2- 
The  subscript  2  represents  conditions  at  the  exit  of  the  oxidizer  duct.  At  the  surface  of  the  liquid  pool  the  temperature  T 1 
can  be  obtained  from  the  equations  describing  the  vapor-liquid  equilibrium.  For  simplicity,  Ti  is  presumed  to  be  equal  to 
the  normal  boiling  point. 

The  strain  rate  characterizing  the  flow  field  is  given  by 


a  = 


(1) 


This  equation  is  obtained  from  an  asymptotic  theory  where  the  Reynolds  number  of  the  laminar  flow  at  the  boundary  is 
presumed  to  be  large  [5]. 

Critical  conditions  of  extinction  are  given  by  the  strain  rate  aF  and  the  oxygen  mass  fraction.  Critical  conditions  of 
autoignition  are  given  by  the  strain  rate  ai  and  the  temperature  at  the  oxidizer  boundary  T2j.  The  extinction  experiments 
are  carried  out  at  room  temperature.  The  fuels  tested  are  n-heptane,  iso-octane,  /(-octane,  n -decane,  /(-dodecane,  n- 
hexadecane,  o-xylene,  JP-10,  diesel,  and  JP-8. 


Experimental  Procedure 

A  detailed  description  of  the  burner  employed  in  the  experiments  is  given  elsewhere  [6].  Computer-regulated  flow 
controllers  measure  the  flow  rates  of  the  gases.  The  calibrated  accuracy  of  these  mass  flow  controllers  is  ±1%.  The 
velocities  of  the  reactants  at  the  boundaries  are  presumed  to  be  equal  to  the  ratio  of  their  volumetric  flow  rates  to  the 
cross-section  area  of  the  duct.  The  temperature  of  the  oxidizer  stream  is  measured  by  a  thermocouple  at  the  exit  of  the 
duct.  The  measured  temperature  is  corrected  for  radiation  by  the  equation  Tg  =  Ttc+(e<Td Ttc4/2A)  A  where  Tg  is  the  gas 
temperature  and  Ttc  is  the  thermocouple  temperature,  a  the  Boltzman  constant  (5.67'lCT8W/m2K4),  A,  the  thermal 
conductivity,  d  the  bead  diameter  of  the  thermocouple,  and  A  the  view  factor. 


Extinction  Experiments 

In  extinction  experiments  the  distance  between  the  oxidizer  duct  and  the  liquid  pool  is  L  =  10  mm.  The  temperature 
of  the  oxidizer  stream  at  the  boundary  was  maintained  at  a  constant  value  of  298  K.  At  given  values  of  oxygen  mass 
fraction  and  a  strain  rate  the  flame  is  stabilized  in  the  mixing  layer  between  the  oxidizer  duct  and  the  surface  of  the  pool. 
The  velocity  of  the  oxidizer  stream  is  increased  until  the  flame  extinguishes. 


Autoignition  Experiments 

In  the  autoignition  experiments  the  distance  between  the  duct  and  the  liquid  pool  is  L  =  12  mm.  The  oxidizer  mole 
fraction  is  held  constant  at  Xo2,2  =  0.21  (mass  fraction  Y>2,2  =  0.233).  The  oxidizer  and  the  evaporating  fuel  are  mixing 
continuously  close  above  of  the  stagnation  plane  (see  Figure  1).  At  a  certain  strain  rate  the  oxidizer  temperature  is 
gradually  increased  until  autoignition  takes  place. 


Results  and  Discussion 
Extinction 

Figure  2  shows  the  experimental  data  as  a  function  of  the  oxygen  mass  fraction  Y>2,2  over  the  strain  rate  at 
extinction  in. 

For  a  given  fuel  type  and  oxygen  mass  fraction,  the  flame  will  extinguish  if  the  strain  rate  of  the  oxidizer  is  higher 
than  eb  E.  That  means  that  for  a  given  value  of  oxygen  mass  fraction  the  region  a2<aE  represents  an  area  of  flammable 
conditions. 
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Figure  2  Experimental  data  showing  the  oxygen  mass  fraction  Y02,2  as  a  function  of  the  strain  rate  aE  at  extinction.  The 
symbols  represent  measurement. 


All  studied  fuels  show  the  same  tendency.  The  strain  rate  at  extinction  increases  with  increase  of  the  oxygen  mass 
fraction.  Also,  the  investigated  fuels  show  a  degressive  behavior.  With  increasing  strain  rate  the  mass  fraction  Yo2i2  also 
increases,  but  the  higher  the  strain  rate,  the  smaller  the  oxygen  mass  fraction  difference  between  two  adjacent 
measurement  point. 

The  order  of  the  reactivity  of  different  fuels  is  nearly  independent  from  the  strain  rate.  Fuels  with  a  smaller  value  of 
C-atoms  are  generally  more  reactive  then  fuels  with  long-chained  carbons.  Not  surprisingly  is  that  adding  a  less  reactive 
fuel  to  another  fuel  decreases  the  overall  reactivity  of  fuel  mixture  as  it  is  done  in  this  study  for  n-decane  mixing  with  o- 
xylene.  With  certain  combinations  it  is  possible  to  achieve  a  fuel  mixture  that  shows  similarly  chemical  characteristics  as 
e.g.  diesel.  That  doesn’t  imply  that  a  found  mixture  with  the  desired  behaviors  for  extinction  shows  the  same  behavior 
for  ignition. 


Autoignition 


Figure  3  shows  the  experimental  data  of  the  oxidizer  temperature  T2ji  as  a  function  of  the  strain  rate  a  at 
autoignition.  These  data  represent  the  critical  conditions  of  autoignition  of  liquid  fuels.  The  symbols  indicate 
measurements. 

At  a  given  value  of  strain  rate,  autoignition  will  take  place  if  the  temperature  of  the  oxidizer  stream  is  greater  than 

T2,i. 

Experimental  data  show  that  the  value  of  T2j  increases  with  increasing  a2.  The  fuel  with  the  highest  autoignition 
temperature  is  o-xylene  followed  by  iso-octane.  The  lowest  oxidizer  temperature  for  autoignition  is  necessary  for 
hexadecane.  That  means  that  hexadecane  is  easier  to  ignite  than  o-xylene  and  applies  also  for  the  order  of  the  fuels  in 
between.  Noteworthy  is  the  crossover  between  n-decane  and  n-dodecane  among  the  other  tested  fuels.  n-Decane  and  n- 
dodecane  ignite  easier  at  lower  strain  rate  than  at  higher  strain  rate  compared  to  the  other  fuel  types.  Also  diesel  and  n- 
octane  show  a  faster  decrease  in  reactivity  over  the  strain  rate  compared  to  the  other  investigated  fuels.  That  means  that 
there  is  an  influence  of  the  strain  rate. 
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Figure  3  Experimental  data  showing  the  oxidizer  temperature  T2I  as  a  function  of  the  strain  rate  a2  at  autoignition.  The 
symbols  represent  measurements. 


Conclusion 

Model  fuels  that  show  similar  chemical  behaviors  as  diesel  can  be  developed  by  mixing  higher  and  lower  reactive 
short -chained  fuels. 

Previous  studies  considered  «-heptane  as  surrogate  diesel.  The  present  study  shows  that  n-heptane  does  not 
reproduce  extinction  and  autoignition  characteristics  of  diesel.  A  mixture  consisting  of  n-decane  and  o-xylene  is  a  quite 
good  approach  to  create  a  model  fuel  for  diesel.  In  both  experiments  the  mixture  of  «-decane  80%  and  o-xylene  20% 
shows  slightly  too  high  values  compared  to  the  results  of  diesel. 
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Appendix  D:  The  influence  of  Water  Vapor  on  the  Critical 
Conditions  of  Extinction  and  Autoignition  of  Premixed 

Ethene  Flames 

M.  Geieregger,  R.  Seiser  and  K.  Seshadri 


1  Introduction 

An  experimental  and  numerical  study  is  carried  out  to  characterize  the  influence  of  water  va¬ 
por  on  critical  conditions  of  extinction  and  autoignition  of  strained  laminar  premixed  flames. 
The  counterflow  configuration  is  employed.  The  fuel  used  is  ethene  (C2H4).  Water  vapor  has 
a  physical  as  well  as  chemical  influence  on  critical  conditions  of  extinction  and  autoignition. 
Addition  of  water  vapor  lowers  the  temperature  in  reaction  zone.  This  makes  the  flame  (pre¬ 
mixed  and  nonpremixed)  easier  to  extinguish  and  more  difficult  to  ignite.  This  is  the  physical 
influence  of  water.  The  chemical  influence  of  water  arises  from  its  chaperon  efficiency  in  three 
body  reactions.  Consider  the  three  body  reaction  H  +  O2  +  M  — ►  HO2  +  M,  where  M  rep¬ 
resents  any  third  body.  This  reaction  is  chain  breaking.  The  concentration  of  the  third  body, 
Cm  appearing  in  the  reaction  rate  of  this  three  body  reaction  is  calculated  using  the  relation 
Cm  =  \pW/(RT)\  Ya=  1  ViYi/Wi  where  p  denotes  the  pressure,  T  the  temperature,  R  the  gas 
constant,  W  is  the  average  molecular  weight  and  Tj,  Wj,  and  rp ,  are  respectively  the  mass 
fraction,  the  molecular  weight,  and  the  chaperon  efficiency  of  species  i.  In  hydrocarbon  air 
flames,  the  value  of  the  chaperon  efficiency  of  water  vapor  is  larger  than  many  other  species 
in  the  reaction  zone.  This  applies  to  many  other  three  body  reactions.  Thus  addition  of 
water  increases  the  reaction  rate  of  three  body  reactions.  A  three  body  reaction  that  is  chain 
breaking  will  make  the  flame  easier  to  extinguish  or  more  difficult  to  ignite  when  water  is 
added  to  the  reactant  mixture. 

Steady,  axisymmetric,  laminar  flow  of  two  counterflowing  streams  toward  a  stagnation 
plane  is  considered.  Figure  1  shows  a  schematic  illustration  of  the  counterflow  configuration 
employed  in  the  present  experimental  and  numerical  study.  The  counterflow  configuration 
comprises  two  ducts.  The  distance  between  the  exits  of  the  ducts  is  L.  Studies  on  premixed 
systems  are  carried  out  by  injecting  a  premixed  reactant  stream  made  up  of  ethene,  oxygen 
and  nitrogen  from  one  duct,  and  an  inert-gas  stream  of  N2  from  the  other  duct.  Water  vapor 
is  added  to  the  premixed  reactant  stream.  The  mass  fraction  of  ethene,  the  mass  fraction  of 
oxygen,  the  mass  fraction  of  water  vapor,  the  temperature,  and  the  component  of  the  flow 
velocity  normal  to  the  stagnation  plane  in  the  premixed  reactant  stream  at  the  exit  of  the 
duct  are  Yy  \,  Yq2,i  ,  Th2o,i,  T\ ,  and  Vj ,  respectively.  The  temperature,  and  the  component 
of  the  flow  velocity  normal  to  the  stagnation  plane  in  the  inert-gas  stream  at  the  exit  of  the 
duct  are  T2,  and  V2,  respectively. 
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Figure  1:  Schematic  illustration  of  the  counterflow  configuration. 


Experimental  studies  are  conducted  with  the  momenta  of  the  counterflowing  streams 
PiVj2,  i  =  1,2  kept  equal  to  each  other.  Here  p  is  the  density.  This  condition  ensures  that 
the  stagnation  plane  formed  by  the  two  streams  is  approximately  in  the  middle  of  the  region 
between  the  two  ducts.  The  tangential  components  of  the  flow  velocities  at  the  exits  of  the 
ducts  are  presumed  to  be  equal  to  zero  (plug-flow  boundary  conditions).  The  value  of  the 
strain  rate,  a,  defined  as  the  normal  gradient  of  the  normal  component  of  the  flow  velocity, 
changes  from  one  duct  exit  to  the  other.  [1],  The  value  of  a  on  the  inert  gas  side  of  the 
stagnation  plane  is  presumed  to  be  given  by  [1] 
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Equation  (1)  is  obtained  from  an  asymptotic  theory  where  the  Reynolds  numbers  of  the 
laminar  flows  evaluated  at  the  exits  of  the  ducts  are  presumed  to  be  large  [1]. 


Description  of  Experimental  and  Numerical  Studies 

A  detailed  description  of  the  burner  is  given  elsewhere  [2-5] .  The  flow  rates  of  gases  are  mea¬ 
sured  by  computer-regulated  mass  flow  controllers.  The  mass  flow  controllers  were  calibrated 
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using  a  wet-test  meter.  The  calibrated  accuracy  is  ±1%.  The  gas  velocities  at  the  exit  of  the 
ducts  are  presumed  to  be  equal  to  the  ratio  of  their  volumetric  flow  rates  to  the  cross-section 
area  of  the  ducts.  The  duct  from  which  the  inert-gas  stream  is  injected,  is  equipped  with 
a  heating  device,  that  allows  the  stream  to  be  preheated.  The  temperature  of  the  heated 
stream  at  the  exit  of  the  duct  is  measured  using  a  Pt-Pt  13%  Rh.  The  measured  temperature 
are  corrected  for  radiative  heat  losses  assuming  spherical  shape  of  the  junction,  a  constant 
Nusselt  number  of  2.0,  and  a  constant  emissivity  of  0.2  [6].  The  accuracy  of  the  corrected 
temperature  is  expected  to  be  better  than  ±25  K.  Water  is  added  to  the  premixed  reactant 
stream  by  flowing  the  air  through  a  heated  bath  of  water.  This  is  done  in  a  vaporizer  similar 
to  that  described  in  Ref.  [7].  After  leaving  the  vaporizer,  the  air  and  water  vapor  are  mixed 
with  ethene  and  nitrogen.  The  amount  of  water  is  calculated  from  the  vaporizer  temperature 
by  assuming  liquid  vapor  equilibrium.  The  performance  of  the  vaporizer  was  tested  with  a 
hygrometer  and  showed  above  99.0%  saturation. 

The  experiments  were  carried  out  at  a  pressure  p  =  1.013  bar.  The  premixed  reactant 
mixture,  made  up  of  ethene,  air,  water  vapor  and  nitrogen,  is  characterized  by  the  equivalence 
ratio,  cj)  =  3Yf,iWo2/0'c>2  jWf)'  The  adiabatic  temperature,  Ta d,  can  be  calculated  using 
the  concentration  of  fuel,  oxygen,  nitrogen  and  water  vapor  in  the  premixed  reactant  stream. 
Experiments  are  carried  out  with  the  concentration  of  various  species  in  the  premixed  reac¬ 
tant  stream  so  chosen  that  the  adiabatic  temperature  is  constant.  This  removes  the  physical 
influence  of  water  vapor  addition  on  critical  conditions  of  extinction  and  autoignition.  This 
way  the  chemical  effect  of  the  addition  of  water  can  be  studied.  In  the  extinction  experiments 
the  temperature  of  the  premixed  reactant  stream,  T\  =  353  K,  and  the  temperature  of  the 
inert-gas  stream,  T2  =  298  K  and  the  equivalence  ratio  is  unity.  The  concentration  of  the 
reactants  are  so  chosen  that  the  adiabatic  temperature  of  the  premixed  reactant  stream  are 
Tj  ad  =  1823  K  and  Tj  a(j  =  1920  K.  In  the  extinction  experiments  a  flame  is  established  in  the 
mixing  layer  between  the  premixed  reactant  stream  and  the  nitrogen  stream.  The  strain  rate 
is  increased  until  extinction  is  observed.  This  is  done  for  mole  fractions  of  water  of  0.0-0.23. 
The  strain  rate  at  extinction  is  o^e-  In  the  autoignition  experiments  the  temperature  of  the 
premixed  reactant  stream,  T\  =  393  K.  The  nitrogen  stream  is  heated  and  the  equivalence 
ratio  is  0.9.  Initially  there  is  no  flame  present  in  the  mixing  layer.  The  nitrogen  stream 
temperature  X2  is  increased  until  autoignition  is  observed.  The  last  value  before  ignition  is 
recorded  as  the  autoignition  temperature  T2.i-  Two  sets  of  autoignition  experiments  are  con¬ 
ducted.  One  with  Tijad  =  1800  K  and  one  with  T\  .A(\  =  1900  K.  The  strain  rate  02  was  held 
constant  at  400  s-1.  The  mole  fractions  of  water  range  from  0-0.18  in  the  2073  K  case  and 
from  0-0.13  in  the  2173  K  case. 

The  numerical  computations  are  carried  out  using  a  computer  program  called  FlameMas- 
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ter,  that  was  developed  at  RWTH- Aachen  [8].  At  the  boundaries  of  the  computational  domain 
the  mass  flux  of  the  reactants  in  the  premixed  reactant  stream,  the  mass  flux  of  nitrogen  in 
the  inert-gas  stream,  and  their  velocities  are  specified  corresponding  to  those  used  in  the 
experiments.  Plug-flow  boundary  conditions  are  employed  in  the  calculations.  The  conserva¬ 
tion  equations  of  mass,  momentum  and  energy  and  the  species  balance  equations  used  in  the 
formulation  of  the  numerical  problem  are  summarized  elsewhere  [8,  9] .  The  species  balance 
equations  include  thermal  diffusion  and  the  energy  conservation  equation  includes  radiative 
heat  losses  from  carbon  dioxide  and  water  vapor  [8].  Buoyancy  is  neglected.  Two  chemical- 
kinetic  mechanism  are  employed  in  the  numerical  calculations.  One  mechanism,  called  the 
San  Diego  Mechanism,  was  developed  at  the  University  of  California  at  San  Diego  [10].  The 
other  mechanism  called  the  LLNL  mechanism  was  developed  at  Lawrence  Livermore  Labora¬ 
tories  [3].  The  San  Diego  mechanism  was  previously  employed  to  predict  various  aspects  of 
premixed  and  nonpremixed  combustion  of  ethyne  (C2H2),  ethene  and  ethane  [11-14].  Chem¬ 
ical  species  containing  three  carbon  atoms  such  as  propene  are  included,  but  considered  as 
intermediate  products.  The  LLNL  mechanism  employed  here  was  developed  to  predict  struc¬ 
ture,  extinction  and  autoignition  of  heptane  flames  [3] .  The  LLNL  was  simplified  by  removing 
reactions  in  which  species  with  three  or  more  carbon  atoms  appear. 

Results  and  Discussion 

Figures  2  and  3  show  critical  conditions  of  extinction  and  autoignition.  The  points  represent 
measurements  and  the  lines  are  results  of  numerical  calculation  obtained  using  the  San  Diego 
mechanism  and  the  LLNL  mechanism.  Figures  2  shows  the  strain  rate  at  extinction,  02, Ej 
as  a  function  of  the  mole  fraction  of  water  vapor  in  the  premixed  reactant  stream,  Xh2o,i- 
Experimental  data  and  results  of  numerical  calculation  with  LLNL  mechanism  show  that  <32, e 
decreases  with  increasing  values  of  Xh2o,i-  Numerical  results  obtained  using  the  LLNL  mech¬ 
anism  agrees  well  with  the  experimental  data.  Numerical  calculations  with  the  San  Diego 
mechanism  show  02, e  to  slightly  increase  with  increasing  values  of  Xh2o,i-  At  a  given  value 
of  Ah2o,1i  the  value  of  02, e  calculated  using  the  San  Diego  Mechanism  is  higher  than  the 
measurements. 

Figure  3  shows  the  temperature  of  the  inert-gas  stream  at  autoignition,  T2.1,  as  a  function 
of  the  mole  fraction  of  water  vapor  in  the  premixed  reactant  stream,  Jh2o,i-  Experimental 
data  and  results  of  numerical  calculation  with  LLNL  mechanism  show  that  T2.1  increases  with 
increasing  values  of  Ah2o,i-  Numerical  calculations  with  the  San  Diego  mechanism  show  TM 
to  slightly  decrease  with  increasing  values  of  -Xh2o,i-  Thus  experimental  data  and  numerical 
calculation  with  the  LLNL  in  Figs.  2  and  3  show  that  the  flame  gets  weaker  with  addition 
of  water,  whereas  numerical  calculations  with  the  San  Diego  mechanism  show  that  the  flame 
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Figure  2:  The  strain  rate  at  extinction,  as  a  function  of  the  mole  fraction  of  water  vapor 
in  the  premixed  reactant  stream,  -Xh2o,i-  The  points  represent  measurements  and  the  lines 
are  results  of  numerical  calculation  using  the  San  Diego  mechanism  and  the  LLNL  mechanism. 
Experiments  and  calculations  are  carried  out  keeping  the  adiabatic  flame  temperature  constant 
at  1823  K  and  1920  K,  respectively. 
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Figure  3:  The  temperature  of  the  inert-gas  stream  at  autoignition,  X^i,  as  a  function  of  the 
mole  fraction  of  water  vapor  in  the  premixed  reactant  stream,  Xh2o,i-  The  points  represent 
measurements  and  the  lines  are  results  of  numerical  calculation  using  the  San  Diego  mechanism 
and  the  LLNL  mechanism.  Experiments  and  calculations  are  carried  out  keeping  the  adiabatic 
flame  temperature  constant  at  2073  K  and  2173  K,  respectively.  The  equivalence  ratio  is  0.9. 
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gets  stronger  with  addition  of  water. 


Calculations  with  the  San  Diego  mechanism  show  that  among  the  three  body  reactions, 
the  elementary  reactions  H  +  O2  +  Mg  =  HO2  +  Mg,  and  HCO  +  M4  =  H  +  CO  +  M4  play 
key  roles  in  influencing  critical  conditions  of  extinction  and  autoignition.  In  the  San  Diego 
mechanism  the  chaperon  efficiency  of  water  vapor  is  Mg  =  7,  and  M4  =  12.0  respectively.  Fig¬ 
ure  4  shows  the  influences  of  these  elementary  reactions  on  critical  conditions  of  extinction. 
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Figure  4:  Numerical  calculated  values  of  the  strain  rate  at  extinction,  o^e,  as  a  function  of  the 
mole  fraction  of  water  vapor  in  the  premixed  reactant  stream,  -Xh2o,i-  The  calculations  were 
performed  using  the  San  Diego  mechanism.  The  figure  shows  the  inffuences  of  the  elementary 
reactions  H  +  O2  +  Mg  =  HO2  +  Mg,  and  HCO  +  M4  =  H  +  CO  +  M4 

Figure  4  shows  that  value  of  02, e  decreases  with  increasing  mole  fraction  of  water  vapor  in 
the  premixed  reactant  stream  if  the  value  of  Mg  is  increased  from  7  to  21,  or  if  the  value  of 
M4  is  decreased  from  12  to  1. 

Figure  5  shows  that  the  values  of  the  temperature  of  the  inert-gas  stream  at  autoignition, 
T2J  increases  with  increasing  mole  fraction  of  water  in  the  the  premixed  reactant  stream  if  the 
Mg  is  increased  from  7  to  21,  or  if  the  value  of  M4  is  decreased  from  12  to  1.  Similar  results  are 
obtained  if  the  value  of  chaperon  efficiency  in  all  three  body  reactions  are  presumed  to  be  unity. 

The  elementary  reaction  H  +  O2  +  Mg  =  HO2  +  Mg  is  chain  breaking.  Therefore  increasing 
the  rate  of  this  reaction  makes  the  flame  easier  to  extinguish  and  harder  to  ignite.  On  the 
other  hand  the  elementary  reaction  HCO  +  M4  =  H  +  CO  +  M4  produces  radicals.  Therefore 
increasing  the  rate  of  this  reaction  makes  the  flame  harder  to  extinguish  and  easier  to  ignite. 
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Figure  5:  Numerical  calculated  values  of  the  temperature  of  the  inert-gas  stream  at  autoigni¬ 
tion,  T2.1,  as  a  function  of  the  mole  fraction  of  water  vapor  in  the  premixed  reactant  stream, 
Jh20.i  •  The  calculations  were  performed  using  the  San  Diego  mechanism.  The  figure  shows 
the  influences  of  the  elementary  reactions  H  +  O2  +  M@  =  HO2  +  M@,  and  HCO  +  M4  =  H 

+  CO  +  m4 


This  is  consistent  with  the  results  of  numerical  calculations  in  Figs.  4  and  5. 
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